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ABSTRACT 

This report describes an interpretive system which transforms the 650 into 
a three-address, floating -decimal, general-purpose computer, primarily suited 
for scientific and engineering calculations. The system is complete in the sense 
that all mathematical, logical and input-output operations normally called for in 
such calculations can be performed within the system, i. e. , without reference to 
the basic operation codes of the 650. The guiding principles in designing the sys- 
tem have been ease of use, as defined in the introduction, high speed of arith- 
metic and frequently used logical operations and full accuracy and range for the 
elementary transcendental functions. 

The report serves a dual purpose. It presents the external characteristics 
of the interpretive system to the potential user by means of detailed explanations 
accompanied by illustrative examples, assuming no previous familiarity with 
internally programmed machines. It also describes the internal structure of the 
system to the professional designer of such systems, enabling him to modify it 
to suit his particular needs or to borrow ideas or building blocks he may find 
useful. 

The system is available in punched card form to anyone who requests it. 
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Note: The material of immediate concern to those who wish to learn how to 
program problems in the interpretive system is contained in sections II-X. Sec- 
tion I is devoted to general considerations and may be bypassed. Section XI deals 
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ested in the design of interpretive systems, but the discussion of possible mod- 
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The experienced programmer may absorb the essentials of the system by 
reading the definitions of the operations. Page references to them are given in 
the summary of operation codes. 
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I INTRODUCTION 



I. 1. GENERAL DESIGN CONSIDERATIONS 



The use of most existing computing devices whose degree of automatic per- 
formance substantially exceeds that of a desk calculator entails certain problems 
not encountered in desk computing. To cope with these problems, one may incor- 
porate additional circuitry into the machine- -this, indeed, appears to be the trend 
in recently announced commercially available machines- -or, alternatively, one 
may program, in terms of the basic language of the machine, a system or super- 
language in terms of which the general user will program his problems. The user 
may consider the machine and the super -language as one entity, and no knowledge 
of the basic machine language is required of him. Before actual calculation, the 
programmer's instructions are translated by the machine into the basic language. 
If this translation or interpretation takes place each time an instruction is to be 
executed, rather than once for all at the beginning of a problem, the super-language 
is referred to as an interpretive language or system. Limitations in storage capac- 
ity may necessitate the choice of an interpretive system rather than a system of 
the once-for-all type in the case of most small or medium- sized computers. 

The designers of an interpretive system are faced with a very large number 
of decisions. To provide a basis of motivation for these decisions, it is convenient 
to list here, in somewhat arbitrary order, some of the above-mentioned problems 
which the present interpretive system proposes to solve. All of them may funda- 
mentally be measured in terms of total time spent by a programmer in learning to 
use the machine and in using it on a specific problem. In this sense, the "ease of 
use" referred to in the abstract above is implicitly defined by the list that follows. 
The price paid for the saving of programmer time is, of course, to be found in 
substantially reduced speed of operation. 

A. Scaling 

The storage medium- -paper- -normally used in desk computing places no 
practical restriction on the size of numbers or on the location of decimal points. 
In using a computer that automatically stores intermediate results in registers 
of fixed length and with the position of the decimal point fixed in advance, a great 
deal of time must in most cases be spent on estimating the range of all intermediate 
results to prevent errors due to overflow. The well-known way of avoiding this 
at the expense of a very substantial increase in the internal complexity of the arith- 
metic operations is to represent each non-zero number in floating -decimal form, 
i.e., as a signed quantity whose absolute value lies in a fixed range, accompanied 
by an exponent of 10 or decimal point indicator. 

B. Length and Complexity of the Program 

Floating decimal arithmetic and frequently needed special functions could 
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be incorporated into a program written in the basic machine language in the form 
of a set of subroutines reached by a two-way transfer of control, ( ,r calling se- 
quence" or "basic linkage") and there are indeed problems for which this is the 
best choice. In many cases—particularly in the case of relatively short prob- 
lems where the results are needed quickly- -a further reduction of the program- 
ming effort is desirable. This may be achieved by combining under single oper- 
ation codes, in an interpretive language, groups of steps in the basic language 
needed for performing frequently occurring tasks. For example, a single instruc- 
tion in which three locations are specified may be used for adding two arbitrarily 
located numbers and storing the result or a block of information of any length 
may be punched on cards as a result of a single instruction. In particular, the 
task of repeating a calculation a specified number of times, each time with appro- 
priate modifications, must be made easy, and the interpretive system described 
in this report goes as far as is believed possible in this direction by providing an 
order {"LOOP") with which simple cases of this task can be handled by a single 
instruction. 

To preserve the simplicity gained by introducing an interpretive system, the 
system must be made complete or self contained so that" most problems can be 
conveniently programmed without reference to two different systems of operation 
codes, although, of course, leaving and re-entering the interpretive system should 
be made possible in order to provide the experienced programmer with complete 
flexibility. 

C. Restrictions 

In desk computing, one cannot fail to notice if the argument for which a func- 
tion value is to be found in a table falls outside the range covered by the table or 
if in transcribing a set of numbers from one area on a piece of paper to another, 
overlapping area, some of the numbers to be transcribed are erased before tran- 
scription. A machine will avoid or detect such blunders only if programmed to do 
so, and as much as possible of this programming should be included in the inter- 
pretive system. For example, it is desirable that all mathematical functions in- 
cluded in the system be available for the full range of argument consistent with their 
definition and with the limitations imposed by the machine itself and that they are 
computed to the full accuracy of the number system used. Error stops indicating 
violations of unavoidable restrictions should be included to the fullest extent that 
space limitations permit. 

D. Program Testing 

The usefulness of a general -purpose computer or interpretive system depends 
decisively on the methods provided for testing ("debugging") new programs --for 
definitions and details, see Sec. VI. In the case of an interpretive system whose 
operating speed*is only one order of magnitude above the speed of card punching, 
testing by means of a tracing routine included in the system compares favorably to 
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console testing, at least in the case of programmers whose familiarity with the 
machine is limited. Either of these methods is thought of as a tool normally used 
only when memory print-outs have been found insufficient. To facilitate testing 
by any of the methods mentioned and to keep programming as concrete as possible, 
the system described in this report assumes that the actual machine location of 
each instruction and stored number is assigned by the programmer. The system 
may, of course, be used in conjunction with regional or symbolic assembly programs, 

I. 2. CHANGES AND ADDITIONS 

Numerous minor changes suggest themselves when the system is viewed in 
the light of the experience gained in designing it; some of them are discussed in 
detail at the end of the report. Corrections of errors not yet revealed must be ex- 
pected. External changes and additions will undoubtedly be proposed after a period 
of use. The present system should, therefore, be considered primarily as a first 
version which each user may consider changing to better suit his needs. Comments 
and suggestions on both internal and external aspects of the system will be greatly 
appreciated. 
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II GENERAL INFORMATION 



II. 1. THE 650 

The IBM 650 is an electronic computer whose basic storage consists of 
a magnetic drum capable of holding 2000 words (numbers) of ten decimal digits 
and sign. The machine is internally programmed, i. e. , the program of instruc- 
tions which the machine is expected to follow is kept on the drum, and the ma — 
chine automatically reads one instruction at a time from the drum, executes it, 
reads another instruction, and so on. Initially, the program is loaded onto the 
drum from punched cards but each instruction is loaded only once, although the 
machine may be expected to execute it many times in the course of a problem. 
Special orders are inserted into the program to cause repetition of prescribed 
sections the desired number of times. In many cases, the instructions exec- 
uted by the machine are changed or modified under control of the program between 
successive executions., This ability to modify its own instructions is one of the 
distinguishing characteristics of an internally programmed machine. 

In the basic model of the 650, all answers are punched by the machine into 
cards, which may be printed on a separate tabulator. 

The 650 is a general -purpose, fixed-decimal machine and any programmer 
may, with the aid of a detailed manual published by IBM, learn to use it as such. 
There are many large problems and many problems of a data processing nature 
for which fixed-decimal operation is definitely indicated, and the programmer is 
asked to give serious consideration to this alternative for all but the very smallest 
problems, since the gain in machine time over floating -decimal operation (ex- 
plained below) may be as high as 10:1. The machine -language programmer may 
relieve himself of many tasks by using the interpretive system for loading, punch- 
ing, calculation of special functions, etc. , provided 1000 storage locations suffice 
for his problem. (See TR OUT, Sec. IV. 1) 

This report describes a system which enables the programmer to use the 
650 as a floating -decimal machine, without being familiar with the fixed-decimal 
mode of operation. Beginning with the next section, all statements will concern 
the system rather than the 650 itself, but it should be borne in mind that if any- 
thing in the system appears restrictive from the viewpoint of a particular appli- 
cation, — storage capacity, speed, card form, word length, etc. --total or partial 
use of basic 650 coding may be the answer. 

II. 2. THE INTERPRETIVE SYSTEM: STORAGE; DATA AND INSTRUCTION FORM 

When the interpretive system is in use 999 ten-digit storage locations, 
numbered 001-999, are unrestrictedly available to the programmer for storing 
instructions and data. The location 000 has a special use ("previous result") 
which will be discussed below. 
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Throughout the system, numbers upon which mathematical operations are 
performed are stored and used in so-called (normalized) floating -decimal form, 
which will be defined as follows: The number zero is written as ten zeros with a 
plus sign ("machine zero' 1 ). Any number A other than zero is expressed as 



A=± A, 



io a i 



where 1<C A X <C 10 and -50^^ 49. In the machine, A is written as the pair 
±(A X , ah where a= a x + 50 and Ai is an eight-digit number with seven decimal 
places. The "machine exponent" a is a two-digit (positive) integer l ocated at the 
right end of the number. Non-zero numbers A not in the range 10~ 50 < I AK10 
cannot be correctly used in the machine, and some of the mathematical operations 
will give an error stop if the result would fall outside this range (see STOPS). Num- 
bers loaded into the machine must also be in the form prescribed above, unless 
special precautions are taken (see Sec. X. 2). 



The system instructions are signed ten-digit numbers of the following 



+ 


°1 


A or 2 


B 


C 



form: 



Here, X is a one-digit operation code and B and C are three-digit addresses. The 
three-digit quantity "A or 2 " is interpreted as an address A if 0^ ^ and as an 
operation code O2 if L = 0. The sign of the instruction is used in connection with 
the LOOP order (see LOOP OPERATIONS). The only difference between the mu- 
tually exclusive 0j_ and 2 operations is that all operations which require three 
addresses have been designated 0^, all others, O2. 

An example will illustrate how the addresses are used in a program: To add 
the number stored in register 200 to the number stored in register 201 and store 
the result in 500, the operation code 0i = 1 (ADD) is used and the instruction reads: 
1 200 201 500. To take the square root of the number in 200 and store the result 
in 500, the 2 -operation 2 = 300 (SQRT) is used: 300 200 500. As will be shown 
in later sections, it is also possible to call out an instruction stored in memory 
and operate upon it, e.g. , increase one of the addresses in it. In storage, no dis- 
tinction is made between instructions and data so that the programmer is free to 
use any memory location for storing an instruction or a number as he sees fit. 

To facilitate explanations, the following notation will be used: The ten-digit 
quantity whose storage location has the address A__will be denoted by A; anal- 
ogously, B will denote the contents of location B. C denotes the result of a calcula- 
tion about to be stored in location C. 

In addition to being stored at C, the result C of any mathematical operation 
(i. e. , arithmetic operations or special functions) and of MOVE 000 and CONS 
(see READING) is automatically stored in the special location 000. If this result 
is needed on the next step, (or, more generally, before it has been replaced by 
the result of a subsequent mathematical operation or MOVE 000 or CONS) calling 
it out by using 000 as an A-address will reduce the execution time in the case of 
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the arithmetic operations. Also, time will be saved in any mathematical operation 
by using 000 as a C-address when C will be needed only on the next step. Execu- 
tion times are discussed in detail in a later section, but it should be emphasized 
that timing considerations only affect the running time of a problem, never the 
correctness of results. All locations are accessible at any time. 

3pecial addresses for obtaining frequently needed numbers, such as, and 
1, are not provided by the system. The programmer should load such numbers 
into locations of his own choosing. 
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Ill MATHEMATICAL OPERATIONS 



III. 1. ARITHMETIC OPERATIONS 

The operations will be introduced in an order chosen to facilitate the learn- 
ing process. Later, a concise summary of all operation codes will be given. 
Alphabetic operation codes are listed in addition to the numerical ones merely to 
facilitate programming; they are not introduced into the 650 and need not be used 
at all. 

The result of each arithmetic operation is rounded. If the result is zero, 
a machine zero is given, i.e., the machine exponent will be 00. An error stop 
occurs if the result of a multiplication or division would fall outside the range of 
the floating -decimal number representation; another error stop detects attempts 
to divide by zero (see STOPS). 

A list of the arithmetic operations follows: 

Numer. Alpha, 

code code Function 



0, = 1 



ADD 



Add (in floating-decimal form) the num- 
ber A stored at A to the number B 
stored at B, store the result C at C 
and 000. Abbreviated: 



A+ B = C 



0i ss 2 



SUB 

MPY 
DIV 

NGMPY 
III. 2. SPECIAL FUNCTIONS 



°1= 3 



X = 4 



X s 5 



Subtract: A - B = C 
Multiply: A . B » C 
Divide : A / B = C 
Multiply negatively: -A . B = C 



The system is intended to give eight-digit accuracy (i. e. , an error less than 
1 in the eighth digit) in computing the special functions included whenever the input 
makes this accuracy possible. For trigonometric functions of an argument exceed- 
ing one revolution and for logarithms of numbers near 1, loss of accuracy follows 
from the mathematical properties of the respective functions and stops (which may 
be bypassed by the setting of a console switch) are provided when this loss exceeds 
two digits. For small values of the argument, where an eight-digit, fixed-decimal 
representation of the sine or arc tangent would contain leading zeros, the floating - 
decimal representation would normally introduce meaningless digits at the right 
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end. To reduce this nuisance to a tolerable level and also make possible trig- 
onometric calculations with extremely small arguments, the formulas sin x = x 
and arc tan x = x are used for }x\< . 0025 and |xj< .001, respectively. Those 
interested will find the methods of computing the special functions described in 
Section XI. 3. 

Aside from the limitations imposed by the above mentioned inherent loss 
of accuracy and by the floating -decimal representation of the result, no restric- 
tions apply to the natural range o f the argument for the special functions. Error 
stops will prevent attempts to take the square root of a negative number or the 
logarithm of a non-positive number. 

The special functions (or, more precisely, elementary transcendental 
functions) arer 

Numer. Alpha. Function 

2 = 300 SQRT TIT = C 

2 ss 301 EXP E e B = C 

2 =s 302 LOG E log e B =• C 

2 =303 SIN R sin B = C, B in radians 

2 = 304 COS R cos B = C, B in radians 

= 305 ART R arc tan B = C, C in radians, lc|<TT/2 

» 350 ABS |BJ = C 

2 = 351 EXP 10 10 B = C 

2 = 352 LOG 10 logio^ = ^ 

= 353 SIN D sin B = C, B in degrees 

2 =: 354 COS D cos B =s C, B in degrees 

2 = 355 ART D arc tan B = C, C in degrees, |c|< 90 

Subdivisions of a degree are expressed decimally, not in minutes and sec- 
onds. 

III. 3. MOVE 000 

In many cases (particularly in connection with the use of subroutines) it may 
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be convenient to be able to call out a number B from B and deposit it in C, as 
well as in 000, without the time-consuming use of a floating -decimal arithmetic 
operation. This is accomplished by the logical operation 0^ 9 ("MOVE") with 
A =r 000. The normal use of MOVE with A £ 000 is described in Sec. IV. 4. 



III. 4. AN EXAMPLE 

For the benefit of anyone with no previous computer experience, a simple 
example illustrating the use of the mathematical operations will be inserted here. 
Suppose that, as a part of a program which is assumed to be already on the drum, 
it is desired to evaluate the function 



f(x) = 



sin x 



/ 



1+ e 



-x3 



Here, x in radians is assumed to be in storage register 500 and the constant 1 in 
600. The quantity e" x3 is to be stored in 501, and f(x) in 502. A program might 
look as follows: 



Alpha. 


°1 


A or O2 


B 


C 


MPY 


3 


500 


500 


000 


NGMPY 


5 


000 


500 


000 


EXP E 





301 


000 


501 


ADD 


1 


000 


600 


000 


SQRT 





300 


000 


400 


SIN R 





303 


500 


000 


DIV 


4 


000 


400 


502 



-X' 



Comments 
x 2 

, store in 501 



1-f. e 



-X" 



J1+ e" x , store 

temporarily 



sin x 



f(x), store in 502 



The extensive use of the "previous result" address, 000, is worth noting, 
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IV LOGICAL, OPERATIONS 



IV. 1. TRANSFER OPERATIONS 

Suppose the machine has been instructed (see LOADING) to begin a program 
by executing the instruction stored in, say, location 101. When this execution is 
completed, the machine will automatically execute instruction 102, then 103, etc., 
until tolcLby the program to do otherwise. Operations whose primary function is 
to influence either the order in which instructions are executed by the machine or 
the selection of stored data upon which the instructions make the machine operate 
will be called logical operations. A simple example of such an operation is 02 = 
203, "Transfer Control". If in the sequence 101, 102, 103 above, instruction 103 
should read t! TR 203 000 080", the next instruction executed by the machine would 
be 080 instead of 104. This may be expressed by saying that "control was trans- 
ferred to 080". The B-address was ignored in this case. The transfer of control 
may be made to depend on the result of calculations (mathematical or logical) in 
which case a "conditional transfer" is said to occur. Logical operations --con- 
ditional or unconditional- -are needed whenever several blocks of instructions, 
located on various parts of the drum, are to be tied together to form a program, 
whenever it is desired to repeat a calculation several times, etc. 

For simplicity in grouping, the following list of transfer operations includes 
two (UNC STOP and NOOP) whose transfer function is of a degenerate nature. In 
a first reading, it may be advantageous to omit the TR SUBR and TR OUT operations. 



Numer. 



2 = 000 



Alpha. 
UNC STOP 



Function 

Stop unconditionally. The machine 
stops regardless of the setting of 
console switches (see CONSOLE) 
and displays 9999 on the address 
lights and B on the display lights. 
This operation should be used only 
where it is intended to discontinue 
the execution of the program, since 
a continuation of the program cannot 
be effected by a simple depression 
of the PROGRAM START key (see 
STOPS). The C-address is ignored 
but should be filled, e. g. , with zeros 
(see LOADING). 



2 = 200 



COND STOP 



Stop conditionally and transfer. The 
machine stops if the PROGRAMMED 
STOP switch on the console is in 
the STOP position. The number 
1120 is displayed on the address 
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Numer. 



Alpha. 



Function 



lights and B on the display lights. 
When the PROGRAM START key 
is depressed, control is trans- 
ferred to C. If the PROGRAMMED 
STOP switch is in the RUN position, 
control is transferred to C with- 
out stopping. 

Caution: If the PROGRAMMED 
STOP switch is on RUN, the stops 
for loss of accuracy in sine, co- 
sine and logarithm and the stop in 
the CONS operation will not occur. 

This operation may be used for stopping at check points in the early running 
stages of a problem, with the option of avoiding the stops during later runs. 



0? = 201 



TR SGN 



2 =: 202 



TR EXP 



Transfer on sign. Control is trans- 
ferred to C if the result of the last 
mathematical operation or MOVE 
000 or CONS is negative, to B if 
it is non-negative, (i. e. , zero is 
regarded as having a plus sign). 

Transfer on exponent. The exponent, 
c, of the result of the last mathe- 
matical operation or MOVE 000 or 
CONS is compared to B (the leading 
digit of B should be 0). Control is 
transferred to C if c > B. If c< B, 
control proceeds to the next instruc- 
tion. 



This operation is particularly suited for the summation of series where terms 
are to be added until they have a prescribed number (50 - B) of leading zeros. For 
example, in order to return to instruction 080 only as long as the absolute value of 
the previous result is . 0001 or greater, one would write "TR EXP 202 046 080". 
This saves a time-consuming floating subtraction preceding the test. The TR EXP 
operation is also intended to take the place of the TR ZERO operation found in most 
systems. Due to the accumulation of small errors during a calculation, it is unwise 
in most cases to expect a result to be exactly zero to eight figures; here a TR EXP 
with a suitably chosen B may prevent a never-ending repetition of a part of a program. 



2 = 203 



TR 



Transfer. Control is transferred to 
C, i.e., the next instruction exe- 
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Numer. 



Alpha. 



2 =204 



TR SUBR 



Function 

cuted will be the one stored at C. 
The B address is ignored but 
should be filled, e. g. , with zeros. 

Transfer to subroutine. The C- 
address of the instruction located 
at C is set equal to B, whereupon 
control is transferred to C. The 
sign of the instruction at C is made 
positive. (For an elucidation and 
applications see SUBROUTINES. ) 



2 =r 205 



TR OUT 



2 = 454 



NOOP 



Transfer out. Control is transferred 
to C and the instruction stored there 
is executed in the basic language of 
the machine (i. e. , outside the inter- 
pretive system). When an instruc- 
tion address 1095 is given in the 
basic language, control is returned 
to the interpretive system beginning 
at the instruction following the TR 
OUT. The B-address of the TR OUT 
instruction is ignored but should be 
filled, e.g. , with zeros. The pro- 
grammer in basic language must be 
careful not to use locations above 999, 
which are occupied by the interpretive 
system. 

No operation. Control proceeds to 
the next instruction. The B- and 
C-addresses are ignored but should 
be filled, e.g., with zeros. 



This operation is likely to occur in connection with tracing (see TRACING, 
particularly the ST TR ERAS operation) or when a superfluous instruction has 
been deleted from a program. 

IV. 2. LOOP OPERATIONS 

A highly repetitive character is required of any problem to be economically 
handled on an automatic computer. In certain instances, such as Newton's iteration 
procedure for the solution of equations, a repetitive process or "loop" is con- 
veniently programmed, merely using conditional transfer operations. In many cases, 
however, some of the instructions to be repeated must be slightly modified in a 



-15- 



systematic way before each new repetition. 



N 
i= 1 



linear expression £ ^ x i with the a^ and the x x 



For example, in the evaluation of a 
stored in blocks of consecutive 



locations, the addresses of a i and xj must be increased by 1 each time a new term 
is to be computed. To facilitate programming of this kind, the system provides 
two methods of so-called address modification. The simpler--but less general-- 
of these methods employs a special counter called the loop box, which is stored 
in a location normally inaccessible to the programmer. If an instruction carries 
a minus sign, the current contents of the loop box will be added to the instruction 
(in fixed-decimal arithmetic and without regard to the sign) before it is executed. 
If, for example, the instruction - 1 531 600 901 is given and the loop box contains 
+ 009 000 009, the instruction actually executed by the machine would read 1 540 
600 910. The original instruction remains unchanged in its storage location. At 
the end of a calculation, an 2 instruction called LOOP enables the programmer 
to increase the contents of the loop box by 1 in one or several address positions and 
to transfer control back to the beginning of the calculation. Hence, the calculation 
may be carried out repeatedly, each time with different addresses used in the exe- 
cution of instructions with minus signs. A test provision included in the LOOP order 
stops the repetition after a specified number of executions and resets the loop box 
to zero for future use. An example will be given after the following list of LOOP 
operations. 



Numer. 
2 = 100 



Alpha. 
LOOP A 



010 



LOOP B 



001 



LOOP C 



Function 

Loop on A. The contents of the 
loop box are increased by 001 
000 000. In other words, the A- 
segment of the loop box is in- 
creased by 1. After the increase, 
the A- segment of the loop box i,s 
compared to the B -address of the 
LOOP instruction. If the A- seg- 
ment is less than B, control is 
transferred to C. If the A- seg- 
ment is equal to B, (or greater, 
which will never be the case in 
normal use) the loop box is reset 
to zero and control proceeds to 
the next instruction. 

Loop on B. Analogous to LOOP A, 
with the B -segment of the loop 
box now being increased and com- 
pared to the B -address of the LOOP 
instruction. 

Loop on C. Analogous to LOOP 
A, with the C- segment of the 
loop box being increased and 
compared to B. 
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Numer. 



2 = 110 



Alpha. 
LOOP AB 



o 2 = 101 



0? =r Oil 



2 = ill 



LOOP AC 
LOOP BC 

LOOP ABC 



Function 

Loop on A and B. Analogous to 
LOOP A. The A- and B- segments 
of the loop box are increased by 
1 and the A-segment is compared 
to B. 

Loop on A and C. Analogous. 

Loop on B and C. Analogous. The 
B -segment is used for the com- 
parison. 

Loop on A, B and C. Analogous. 
The A-segment is used for the 
comparison. 



To illustrate the use of a LOOP order, consider the evaluation of the linear 
expression L(x) = i£ a- x-, where the a. and the x- are stored in memory. In 

i c 1 1 

choosing storage locations for numbers, it is wise to plan in advance how they are 
to be used in the program. In this case, since the a^ and the x. are to be reached 
using the LOOP operation, it is advantageous to store them in blocks of consecutive 
locations, say the a i in 800+ i and the x i in 900+ i, (i - 1, 2, . . . , 20). Suppose 
L(x) is to be stored in 700. For simplicity, assume that register 700 contains zero 
at the beginning of the calculation and that the loop box has been reset. The entire 
program for this calculation might be written as follows: 



Instr. 


No. 


101 




102 




103 




104 





Alpha. 



Sign. 



0. 



A or 0. 



B 



801 


901 


000 


000 


700 


700 


110 


020 


101 



MPY — 3 

ADD 4- 1 

LOOP AB ■+ 

Next instruction in the problem. 

Note that the B -address of the LOOP order simply indicates the number of times 
the arithmetic calculation is to be performed, including the first time when the 
addresses are actually unmodified (modified by adding zero). The practice of 
starting the instruction numbering at, e.g. , 101, rather than 001 facilitates later 
additions to the beginning of a program. 

The loop box is automatically reset at the beginning of a new problem (see 
LOADING), and whenever a transfer out of a loop is effected by a loop order (as 
stated in the definitions above). Hence, the resetting of the loop box need not 
concern the programmer under normal conditions. If the need for resetting the 
loop box should arise, however, this is easily done by giving, e.g., the order 
LOOP A with the B -address 000. According to the definition of LOOP A, this will 
cause control to proceed to the next instruction with a resetting of the loop box. 
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The C-address is irrelevant in this case. This situation would arise if control 
were transferred out of a loop in the middle of it by one of the conditional transfer 
operations. 

It is worth observing that a LOOP operation may be advantageously used in 
some cases where address modification is not involved, simply to repeat a sequence 
of steps a prescribed number of times, e. g. , each time adding a fixed increment to 
a parameter. In such a case, any one of the loop orders could be chosen, (see 
EXECUTION TIMES, however) and no negative instructions would occur. 

The advantages of the loop-box method are its simplicity and high speed and 
the fact that the original instructions remain unchanged in memory. It is limited 
by the fact that there is only one loop box and hence, all instructions to be modified 
are modified in the same way. To handle situations more complicated than this, 
the system provides a set of operations described in the next section. 

IV. 3. ADDRESS CHANGE OPERATIONS 

Many problems can be completely programmed without the use of address 
change operations, and for someone approaching the field of internal programming 
for the first time, it might be advantageous to ignore this section until the need 
for more general logical operations arises. 

The functions of the address change operations are: (a) To increase or 
decrease a designated address of an instruction in storage by any given amount; 
(b) To set such an address to a given value (without reference to its previous value); 
and (c) To transfer control as a result of comparing such an address to a given num- 
ber. 

There are nine (^-operations among the address change_ operations. In each 
of these, the B -address gives the location of the instruction (B) to be changed and 
the C-address is the amount of change. For example, suppose the instruction 
600 750 005 (using the operation 2 = 600,. ADD A) is given and suppose location 
750 contains the instruction 1 320 400 000. Then the A-address, 320, of this in- 
struction will be increased by 005 and the resulting instruction 1 325 400 000 stored 
back in 750. Similarly, if 050 750 333 were given, (using 2 = 050, SET B) the 
instruction in 750 would be changed to read I 320 333 000. In brief: 



Numer. 



2 s 500 



Alpha. 
SET A 



Function 

Set the A-address. _ The A-address 
of the instruction (B) specified 
by Bis set equal to C. 



2 = 050 



SET B 



Set the B~address._ The B-address 
of the instruction (B) specified by 
Bis set equal to C. 



2 s= 005 



SET C 



Set the C-address. The C-address 
of the instruction (B) specified by 
Bis set equal to C. 



-1 ft. 



Numer. 



2 = 600 



Alpha. 
ADD A 



Function 

Add to the A-address. The A- 
address of the instruction (B) 
specified by B is increased by 
C. 



0„ = 060 
2 



ADD B 



Add to the B -address. The B- 
address of the instruction (B) 
specified by B is increased by 
C. 



2 = 006 



ADD C 



Add to the C-address. The C- 
address of the instruction (B) 
specified by B is increased by 
C. 



2 = 700 



o 2 = 



SUB A 

v 



2 = 070 



007 



SUB B 



SUB C 



Subtract from the A-address. 
The A-address of the instruction 
(B) specified by B is decreased 
by C. 

Subtract from the B -address. 
Analogous to SUB A. 

Subtract from the C-address. 
Analogous to SUB A. 



The sign of the instruction being modified remains unchanged and does not affect 
the outcome of the modification. Attempts to increase an address beyond 999 or 
decrease it below will result in erroneous operation not prevented by error stops, 

Three 0^ -operations, TR A, TR B and TR C, complete the set of address 
change operations. In each of them, the A-address specifies the instruction (A) 
to be called out and the B-address is the constant to which a specified address is 
to be compared. In case of inequality, control is transferred to C. For example, 
if the instruction 6 750 325 200 (using 0j = 6, TR A) is given, control will be trans- 
ferred to 200 if the instruction in 750 reads 1 320 400 000 but control will proceed 
ahead if 750 contains 1 325 400 000. Summarizing: 



0, = 6 



TR A 



Transfer on the A-address. The 
A-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 



19- 



Numer, 



°1 = 7 



0, = 



Alpha. Function 

TR B Transfer on the B-address. The 

B -address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 

TR C Transfer on the C-address. The 

C-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 

As an introductory example, the summation in the section on LOOP OPERA- 
TIONS will be programmed again using address change methods. This would be 
an inefficient choice in an actual problem, but it will best illustrate the difference, 
as well as the analogy between the two methods. It is again assumed that register 
700 contains zero at the start, but the steps analogous to the resetting of the loop 
box will be included. 



Inst. 



Alpha. 



Sign 



0. 



A or 0, 



B 



101 


SET A 


-h 





500 


103 


801 


102 


SET B 


+ 





050 


103 


901 


103 


MPY 


+ 


3 


[ ] 


[ ] 


000 


104 


ADD 


+ 


1 


000 


700 


700 


105 


ADD A 


+ 





600 


103 


001 


106 


ADD B 


+ 





060 


103 


001 


107 


TR A 


+ 


6 


103 


821 


103 


108 


Next instruction 


in the i 


problem. 









The brackets in the A- and B -addresses of instruction 103 are used to indicate that 
these addresses are variable and will be supplied by the program before the instruc- 
tion is executed, hence what is written there when the program is loaded into the 
machine is irrelevant. At the end of the program when instruction 108 is reached, 
memory location 103 will contain + 3 821 921 000. It is assumed that the summation 
just programmed is part of a larger problem in which it is used repeatedly. This 
is the reason for the SET A and SET B instructions. If 801 and 901 were simply 
loaded into their respective positions in instruction 103 initially, the summation 
would be performed correctly the first time it is used, but the next time when the 
summation is called for, instruction 103 would read + 3 821 921 000 and erroneous 
calculations would result. The SET instructions could, of course, have been in- 
serted after the completion of the summation, restoring instruction 103 to its proper 
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value for the next application, but this procedure is not recommended because it 
makes it more difficult to restart the problem from the beginning without reload- 
ing the program in case of an interruption (e. g. , error stop) during the loop. 



A more realistic example of the use of address change methods would be a 
calculation involving more than one summation index or parameter. Then, one 
of the fast and convenient LOOP orders would normally be used in the "inner loop", 
i.e., the loop occurring most frequently, with address change operations con- 
trolling the "outer loop" or loops. Suppose, for example, that it is desired to 
calculate S^ ~ 10 



'J 



_ a..x- for j a 1, 2, ... 5, where the a., are stored in 800 + 10. + i 
.** J 1 1 J 1 j 

i=*l 

. ; a ? . in 821, a ? in 822, and so forth), the 

DredT in 700-4- j. It will be assumed that register 

500 contains zero. For completeness, the setting of all variable addresses to their 

initial values for repeated use of the summation program will be included. 



(i. e. , a^j is in 811, a, 2 * n 812, etc. , «. 2 - -« , - ? , 

x^ in 900 + i, and the S. are to be stor« 
1 J 



Instr. 


Alpha. 


Sign 


101 


SET A 


4- 


102 


SET C 


+ 


103 


MOVE 


+ 


104 


MPY 




105 


ADD 


+ 


106 


LOOP AB 


+ 



107 



110 



111 



MOVE 



+ 



TR C 



+ 



A or 0, 
500 
005 
000 

C 1 

000 
110 

000 



107 



Next instruction in the problem. 



400 



108 


ADD A 


+ 





600 


104 


109 


ADD C 


+ 





006 


107 



706 



B 


C 


Comments 


104 


81 r 


Set variable 
' addresses to their 


107 


701 
J 


initial values 


500 


400 


Set register 400 
to zero 


901 


OOO* 1 


"Inner loop" 


400 


400 ' 


r i.e., summation 


010 


104 


on i 



c 1 



01 ol 
001J 

103"! 



Store the result 

Increase addresses 
for next repetition 
in the outer loop 
(j-loop) 

Test for end of 
j-loop 



A superficial examination of this program might suggest that only 1/5 of the 
program is devoted to actual arithmetic calculation (!), but it should be observed 
that in terms of the number of instructions executed by the machine when one 
complete summation is performed, the arithmetic ones are still in the majority, 
and in terms of execution time they comprise about 3/4 of the program. 
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In programming problems involving several loops, it may be helpful to 
consider the structure of a loop in terms of four phases: 

1. Initialization. Where addresses in the loop are set to their ini- 
tial values, registers used for summation are set to zero, etc. . 
The automatic resetting of the loop box and the fact that ad- 
dresses remain unchanged in memory tend to reduce the ini- 
tialization when the loop is controlled by a LOOP operation. 

In the summation problem above, steps 101 and 102 constitute 
the initialization for the outer loop, step 103 is the initialization 
for the inner loop. Notice that step 103 is repeated as a part of 
the outer loop. 

2. Execution. Comprising the mathematical operations of the loop, 
as well as any logical operations associated with a loop inside 
the one being executed. Above, the execution of the inner loop 
consists of steps 104 and 105 and the execution of the outer 
loop consists of 103-107. 

3. Modification. Where addresses, parameter values, etc. , are in- 
creased or decreased. The modification of the inner loop above 
is included in the LOOP instruction. The modification of the 
outer loop consists of steps 108 and 109. The position of the 
modification in the program in relation to the execution and test 
is frequently subject to choice. 

4. Test. Determining whether the loop is completed or further rep- 
etition is required. The LOOP instruction includes the test for 
the inner loop and step 110 is the test for the outer loop. 

Note: It is important to write loops in such a way that all ini- 
tializations are performed by the program, not by loading. If 
this rule is not followed, it will not be feasible to restart the pro- 
gram during testing or after a machine stop without reloading. 
For example, if a register is used for summation, it should be 
reset before being used in the summation loop by moving zero into 
it from another location, not by loading zero into it from a card. 

Many programmers find it helpful in programming a large problem to draw 
a block diagram or flow chart with one box representing each phase of each loop 
and arrows connecting the boxes showing the flow of control. 

The address change operations, particularly the SET operations, are fre- 
quently useful in non- repetitive situations as well. An example of this will be 
found in the section on SUBROUTINES. 

If a program appears to require a large amount of address modification 
and particularly, if this occurs because a quantity whose address is subject to 
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change is needed in many places in the execution of a loop, it may be advantageous 
to write the execution largely in terms of fixed addresses and perform the mod- 
ifications by moving data. Instruction 107, in the example above, illustrates this 
in a simple way: If the registers 700+ j themselves had been used in the summation 
process, (step 105) both the B~ and C-addresses of step 105 would have required 
modification in the outer loop, as well as the C-address of instruction 103. For 
cases where several numbers are to be moved at the same time, a more general 
MOVE operation than the MOVE 000 used so far is available and will be described 
in the next section. 

IV. 4. MOVE 

The MOVE operation is defined as follows: 



Numer. 



°i = 



Alpha. 



MOVE 



Function 

Move. If A £ 000, the block of 
A consecutive words beginning at 
B is moved into the set of A con- 
secutive locations beginning at 
C. The words in the original 
locations are not destroyed, 
except where the two regions 
overlap. The number in location 
000 ("previous result") is not 
affected when C ^ 000. Both 
C y B and C ^ B are permissible. 
An error stop occurs if C 4- A — 
1_^1000. If A= 000, the word 
(B) specified by the B -address is 
moved into location C and into 000, 
It also remains in location B. 



MOVE with A = 000 differs from MOVE with A= 001 only in that the execution time 
with A= 000 is shorter and that the previous result location is affected. 

Note: If a number is to be moved from location B into 000 for use in a TR 
SGN or TR EXP operation on the next step, MOVE 9 001 B 000 must not be used, 
since these transfer operations work strictly according to their definitions (see 
Sec. IV. 1. ). The correct instruction would be MOVE 9 000 B 000. (Internally, 
these transfer operations inspect a duplicate "previous result" location rather than 
000!) 



Xc in 7 05 and the instruc- 



tion MOVE 9 005 701 703 is given. Then x x will be found in 703, x in 704, 



As an example, suppose xj is in 701, x^ in 702, 
[OVE 9 005 701 703 i 
X5 in 707, after execution. 

In conclusion, it should be pointed out that the use of the logical operations 
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is by no means restricted to the straightforward functions for which they are 
primarily intended. The programmer will find innumerable ways of increasing 
the efficiency and elegance of his programs by unusual applications, particularly 
of the address change operations. As a weird example, suppose it is desired 
to multiply the numbers located in registers 1, 4, 9, 16, 25, 36, . . . , 400 (!) 
by a constant located in 600 and store the results in 501, 502, 503, . . . , 520: 



instr. 


Alpha. 


Sign 


°1 


A or 2 


B 


C 


898 


SET A 


+ 





500 


900 


001 


899 


SET C 


+ 





005 


901 


003 


900 


MPY 


— 


3 


[ ] 


600 


501 


901 


ADI? A 


+ 





600 


900 


C ] 


902 


ADD C 


+ 





006 


901 


002 


903 


LOOP C 


+ 





001 


020 


900 


904 


Next instruction in the ■ 


problem. 
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V INPUT-OUTPUT OPERATIONS 



V. 1. CARD FORM 

By a card form is meant a specific assignment of card columns to form fields 
for data, instructions, identification, etc. , in connection with a given program or 
interpretive system. In the 650, information is transmitted to and from cards 
through a control panel, and anyone whose needs call for a special card form can 
adapt it for use in connection with the interpretive system merely by simple control- 
panel wiring. For most needs, the following card form, associated with the inter- 
pretive system, is likely to be found adequate. At this point, only brief definitions 
of the card fields will be given for reference in subsequent sections where their 
use will be explained in detail: 



Columns 



Definition 



1-4 
5-6 
7-9 

10 

11 
12-21 

22 
23-32 

33 
34-43 

44 
45-54 

55 
56-65 

66 
67-76 
77-79 

80 



Card number 

Deck number 

Location 

Word count 

Sign of word 1 

Word 1 

Sign of word 2 

Word 2 

Sign of word 3 

Word 3 

Sign of word 4 

Word 4 

Sign of word 5 

Word 5 

Sign of word 6 

Word 6 

Problem number 

Tracing identification 



The same card form is used in all input-output operations, as well as in trac- 
ing. Both instructions and data are signed ten-digit words and are entirely indistin- 
guishable in connection with input -output operations. 

V. 2. PUNCHING 

At any point in the problem, the machine may be ordered to punch into cards 
the contents of any set of memory locations, together with appropriate identification. 
In some problems, it may be desirable to punch out answers one at a time, perhaps 
together with the values of relevant parameters; in others it may be preferable to 
punch out a large amount of information at less frequent intervals. There are also 
cases where it is advantageous to punch out instructions: In connection with testing 
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(see PROGRAM TESTING) in order to examine a program interrupted at a chosen 
point, and in connection with loading, (see LOADING) in order to reduce the size of 
a deck of cards. All of these ends are served by the following instruction: 



02=^410 PCH Punch cards. The block of consecutive 

words beginning at B arid ending at C 
(inclusive) is punched into cards. Five 
words and a word count of 5 are punched 
into each card but the last, whose word 
count will be the remainder when C-B-4* 
1 is divided by 5. On each card, the 
location from which word 1 was punched 
is punched into columns 7-9. The 
words in storage are not destroyed. A 
cumulative count of the number of cards 
punched during the problem (i. e. , since 
LOADING) is punched into columns 1-4. 
The problem number (see LOADING) is 
punched into columns 77-79 and zero is 
punched into columns 6 and 80. An error 
stop occurs if B ^ C. 

If it is desired to punch six words to a card, this may be done by adding a 
special card behind the punching deck (see LOADING). This card should have an 
x-punch in column 5, 1969 in columns 6-9» 1 in column 10, a 12 punch in column 
11, and 00 0006 0000 in column 12-21. 

The punched cards are likely to be used for one (or both) of two purposes (in 
addition to possible processing on other equipment): The information on thern may 
be printed on a tabulator or they may be loaded (or READ) into the 650 at a later 
time. Details of the printing will not be given here, since they depend on the char- 
acteristics of the tabulator, but the printing form may be assumed to be roughly 
identical with the card form with proper spacing between words. (Suggestions on 
tabulator wiring are given in Sec. XI. 5. ) It is assumed that the suppression of the 
superfluous words punched into the last card, if its word count is not 5, will be 
performed on the tabulator control panel. If this is not feasible, it may be done in 
the 650 by adding three cards to the punching deck. For details, see Sec. XI. 1. 

Selective spacing between lines in printing may be accomplished in several 
ways, even though no operation in the 650 is provided for this purpose. A brief 
discussion will be given here, _since spacing considerations may affect the use of 
the PCH operation in programming. -Through the setting of switches on the tab- 
ulator, a choice of any of the following spacing alternatives maybe provided: 

(a) Single or double spacing. 
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(b) Spacing between every n lines (with n chosen by wiring, nor- 
mally, e. g. , n = 10). 

(c) Spacing after any line whose word count is less than the word 
count of the preceding line. 

(d) Spacing before any line whose location number has a units digit 
smaller than the units digit of the location number of the pre- 
ceding line. 

Alternative (c) is suited for the printing of information punched from fairly 
large blocks of locations by one PCH order. Spacing will occur after each block, 
unless the block length is a multiple of 5, which can be avoided by programming. 
Alternative (d) is intended for information punched repeatedly from the same set 
of locations and provides the option of spacing when the loop is interrupted, e. g. , 
for changing a parameter value, 

V. 3. LOADING 

When a program has been written, and careful inspection reveals no further 
errors, it is key punched into cards following the card form given in Sec. V. 1. To 
reduce to a minimum the number of errors to be found with the aid of the 650, the 
cards should be run through a verifier operated by another person or, alternatively, 
key punched independently by two operators and compared on a reproducer. The 
programmer has the option of specifying the number of words to be punched to a 
card: Punching 5 or 6 to a card will keep the program deck small from the outset 
and eliminate the need for condensing the deck on the 650 later. Punching one word 
to a card is felt by some programmers to facilitate changes. Each card must have 
in columns 7-9 the location into which word 1 is to be loaded, and in column 10 the 
number of words to be loaded from the card into consecutive locations. Columns 
1-6 and 77-80 are not read by the 650 (except that the problem number is read from 
the last card, see below) and may be used by the programmer as he deems best. 
Each column of each field to be used by the machine must contain one and only one 
punch and an error stop is provided to enforce this rule. A 12-punch is used for 
plus, an 11- or x-punch for minus and a 0-punch--not a blank column--for zero. If 
the word count is less than 6, unused word fields and sign columns may be left 
blank. No distinction is made between data and instructions in key punching and load- 
ing. 

LOADING is the process of feeding data and instructions into the machine at 
the beginning of a problem. If the previous user of the 650 was not using the inter- 
pretive system or if there is any reason to doubt that the system is correctly 
stored on the drum, the program deck should be preceded in loading by a deck which 
loads the interpretive system (in 51. 9 seconds) into the memory locations above 999. 
Before the program deck, the programmer may also place a Reset Memory Card, 
which will (in 6. 3 seconds) reset each of the memory locations 001-999 to minus 
zero. (This is useful in connection with the punching out of sections of memory in 
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testing. ) Immediately behind the program deck--no blank cards are used in the 
card reader in connection with this interpretive system- -the programmer places 
one of two nine -card decks to inform the machine whether he wants normal opera- 
tion or TRACING described in a later section. (If he knows that he wants the 
same mode of operation as the previous user, he can omit these cards but the 
gain is only 2. 7 seconds. ) Last, he must place a so-called transfer card with a 
zero punched in column 10, the problem number in 77-79 and the location of the 
instruction at which the program begins in columns 7-9. The word fields on this 
card may be left blank. 

The loading program automatically resets the loop box, the card counter 
(see PUNCHING) and location 000 to zero. 

The order in which the program cards are loaded is irrelevant, unless the 
same location is loaded into from more than one card, in which case the last such 
card, of course, determines the contents of the location. This may occur in 
connection with changes of a temporary nature, which may be placed at the end of 
the deck and later removed, leaving the program in its original form. In the deck 
which loads the interpretive system, the order of the cards must be preserved, 
and an error stop is provided to insure this, thereby ascertaining that no part of 
the system is missing. 

In summary, a complete deck to be loaded must contain: 

System deck (173 cards) 
Reset Memory card (optional) 
Program deck 

Mode-of-operation deck (9 cards) 
Transfer card 

The control console of the 650 need be of almost no concern to the user of 
the interpretive system under normal conditions. He must only make sure that 
all switches on the console are set in a fixed manner required by the system, and 
these settings will now be listed without any description of the function of the 
switches. Certain ways of using the console are described in the sections on READ- 
ING and PROGRAM TESTING. 

Switches Settings 

Storage entry 70 1951 1333 +■ 

Programmed stop Stop (see COND STOP) 

Half cycle Run 

Address selection 1338 (see STOPS) 

Control Run 

Display Upper Accumulator 

Overflow Stop 

Error Stop 

To start a problem, the deck to be loaded is placed in the card reader, and 
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the following keys are depressed in order: 

(1) COMPUTER RESET 

(2) PROGRAM START 

(3) START 



(on the console) 
(on the console) 
(on the card reader) 



When the last card leaves the hopper, the machine stops and the key labelled 

(4) END OF FILE (on the card reader) 

is depressed. If the deck has been correctly put together, the execution of the pro- 
gram will then start automatically. 

The program deck may be run out at any time after loading by depressing the 
START key, unless a READ instruction is contained in the program. Blank cards 
should be inserted into the PUNCH hopper and the START key on the punch side 
depressed. 

To make the 650 produce a condensed program deck in case the program was 
originally key punched one instruction to a card, a PCH instruction should be given 
at the very beginning of the program. This instruction may be bypassed during 
subsequent executions of the program merely by changing the location number on 
the transfer card. 



V.4. READING 

In some problems, particularly in applications of a data processing nature, 
it may be desirable to read information into the machine during the execution of 
the program without manual interference. This is accomplished by the READ 
operation: 



2 = 400 



READ 



Read cards. The block of consec- 
utive storage locations beginning 
at B and ending at C (inclusive) 
is read into from cards. The ad- 
dress B must appear in the location 
field on the first card, as well as 
in the READ instruction, and the 
location field on each card follow- 
ing must contain the sum of the 
word count and location on the 
previous card. The sum of the 
word counts of all cards to be 
read must be C-B-f 1. Violations 
of these requirements, which have 
been included for the programmer's 
protection, will result in error 
stops. 
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The cards to be read should be placed in the hopper of the card reader imme- 
diately following the transfer card (no blank cards). 

The decisions made with the aid of conditional transfers and other logical 
operations are normally based on criteria predetermined by the programmer and 
incorporated into the program. If the programmer wishes to influence the pro- 
gram during its execution, e. g. , on the basis of a result displayed on the console 
in connection with a COND STOP instruction, he may do so using the CONS opera- 
tion: 

2 — 401 CONS Read console. The machine stops 

if the PROGRAMMED STOP switch 
is on STOP. Zero is displayed on 
the display lights and 1131 on the 
address lights. When the PROGRAM 
START key is depressed, the num- 
ber entered on the STORAGE 
ENTRY SWITCHES is stored in 
location C and in 000 (the "previous 
result ,r location). If the PRO- 
GRAMMED STOP switch is on 
RUN, the storing takes place with- 
out a stop preceding. The B-ad- 
dress is ignored but should be 
filled, e.g., with, zeros. 

An example of an application of CONS might be the feeding in of an "educated 
guess" for a starting value in connection with the solution of algebraic equations. 
Another application, involving only the storage entry SIGN switch, might be to 
continue a program until another user is ready to take the machine, at which time 
a change in the SIGN switch setting, interpreted by a TRSGN operation, causes the 
program to punch out intermediate results for later restart. 
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VI PROGRAM TESTING 



VI. 1. MEMORY PRINT-OUT 

The choice of methods for testing ("debugging") a program by comparing 
results of machine calculation to known quantities or to results of independent 
calculations by other means is governed by the relative availability of machine 
time and programmer time. If machine time is freely available, testing with the 
aid of the control console is highly efficient, as well as instructive and enjoyable, 
as soon as a certain facility for operating the console has been acquired. Par- 
ticularly in the case of small problems, the method of tracing --where a card is 
punched for each instruction executed, showing all numerical and logical quantities 
associated with the execution- -may be the most desirable in that it gives an almost 
certain clue to the difficulty within a predictable, if not very short, period of 
machine time and allows the programmer to study the material at his leisure. 
The method most economical of machine time and yet frequently sufficiently 
illuminating is that of memory print-out. It might be suggested that on most prob- 
blems in a busy .but not heavily over-loaded installation, the methods be used in 
the order reverse to that in which they were mentioned here. Some directions 
for their use will now be given. 

The memory print-out method simply consists of inserting temporarily into 
the program at one or several suitably chosen points PCH orders (see PUNCHING) 
calling for the punching of blocks of information- -data or instructions --which, when 
printed on the tabulator, will give a picture of the progress of the program. Since 
1000 words may be punched 5 to a card in two minutes, it is not out of the question 
to punch out the contents of every register used in a problem- -including all the 
instructions --several times. To get the most benefit from this method, the pro- 
grammer should, in any problem that does not threaten to fill the entire available 
memory, avoid using the same storage location for storing different quantities at 
different times whenever feasible, so that as many partial results as possible are 
preserved for the memory print-outs. Whenever a test case of a problem is run, 
even if memory print-out is not chosen as the primary testing method, it would 
certainly be advisable to make the last instruction of the test deck punch out the 
entire memory used. A flexible alternative would be to have scattered through the 
program CONS--TR SGN combinations which transfer control to a PCH order if 
the storage entry sign switch is turned to minus. 

Temporary instructions may be inserted into a program in two ways: Either 
they are included in the normal sequence of instructions when the program is 
initially written and replaced either by NOOP instructions (see Sec. IV. 1) or by 
transfer to the next non-temporary instruction when no longer needed,or else one 
of the regular instructions of the program is replaced by a TR to a vacant location 
L, the regular instruction is placed in L, the temporary ones inLf 1, L + 2, etc. , 
and at the end of this temporary sequence a TR back to the normal program is given. 
In either case, the temporary instructions may (as suggested in LOADING) be kept 
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as a separate deck at the end of the program deck, eliminating any changes in the 
main program deck and simplifying bookkeeping. 

VI. 2. TRACING 

If the tracing deck of nine cards is loaded with the program deck, (see LOAD- 
ING) the machine will automatically start tracing from the beginning of the program, 
as specified by the transfer card. Before the execution of each instruction, a card 
with the following information will be punched: 

Columns Definition 

1-4 Card number (cumulative) 

6 Zero 

7-9 Location of the instruction 

about to be executed. 
10 Six 

11-21 The instruction as stored in 

memory. 
22-32 The instruction as modified 

for execution (i. e. , with the 
loop box added if minus). 
33-43 The contents of the loop box. 

44-54 A if A 4 000, zero if A= 000. 

55-65 B 

66-76 The contents of location 000 

(i. e. , the result of the last 
mathematical MOVE 000 or 
CONS operation). 
77-79 Problem number. 

80 Eight (used by the tabulator 

for automatic selection of a 
different printing form for 
trace cards). 

The punching rate will be 100 cards per minute except in the case of very time-con- 
suming operations, such as, the moving of a large block of information. The ad- 
vantage of punching the trace card before execution is that information will be 
punched for an instruction whose execution is interrupted by an error stop. In the 
case of instructions (such as LOOP or TR EXP) whose B -address does not refer 
to a memory location, the quantity B is irrelevant. Tabulator wiring to suppress 
the printing of B in such cases can be provided if sufficient selection equipment is 
available. The PCH operation is bypassed when the machine is operating in the 
tracing mode, i. e. , PCH is equivalent to NOOP. 

If a program is too long to be traced in its entirety or if this is unnecessary, 
selective tracing may be effected by using the following operations: 
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2 = 450 



START TR 



°2 = 451 



STOP TR 



Start tracing. If the nine-card 
tracing deck has been loaded, 
the machine will start tracing 
from the next instruction. ..If it 
is already tracing, it will con- 
tinue to trace. The B and C 
addresses are ignored. If the 
deck for normal operation has 
been loaded, START TR will 
be equivalent to NOOP. 

Stop tracing. If the machine is 
tracing, it will discontinue trac- 
ing immediately. If it is not 
tracing, STOP TR will be equiv- 
alent to NOOP. The B and C 
addresses are ignored. 



° 2 = 



452 



ST TR ERAS 



Start tracing and erase itself. 
If the tracing deck has been 
loaded, the machine will start 
tracing from the next instruc- 
tion. If it is already tracing, 
it will continue to trace. If the 
deck for normal operation has 
been loaded, . tracing will not 
begin. In all cases, the ST TR 
ERAS instruction will be replaced 
in memory by a NOOP (0 2 = 454) 
during its first (and only!) execu- 
tion. The B and C addresses are 
ignored. 

The bypassing of the PCH operation is in effect as long as the trace program is on 
the drum and is not affected by the selective tracing orders. To make PCH opera- 
tive, the nine-card deck for normal operation must be loaded. 

The purpose of the ST TR ERAS operation is to make it possible to trace the 
repetitive steps of a loop either once or twice and then stop tracing until the loop 
is completed. To get the steps traced once, one may place the pair STOP TR, ST 
TR ERAS at the beginning of the repeated portion of the loop; to get them traced 
twice, one places this pair of instructions at the end immediately preceding the 
test. As a specific example, suppose it is required to trace twice the steps of the 
loop programmed in the section on LOOP OPERATIONS and suppose vacant loca- 
tions are available from 900 up. Assume that the machine is tracing as it enters 
the loop. The original program reads as follows: 
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801 


901 


000 


000 


700 


700 


110 


020 


101 



101 MPY — 3 

102 ADD -+- 1 

103 LOOP AB .'+ 

104 Next instruction in the problem. 

The following instructions could be added as a temporary deck at the end of the pro- 
gram deck: 



102 


TR 


+ 





203 


000 


900 


900 


ADD 


+ 


1 


000 


700 


700 


901 


STOP TR 


+ 





451 


000 


000 


902 


ST TR ERAS 


+ 





452 


000 


000 


903 


TR 


+ 





203 


000 


103 



Notice that the TR instruction gets loaded into 102 after the regular program, replac- 
ing the ADD instruction, as explained in LOADING* This example is, of course, 
unrealistic in that selective tracing would hardly be needed for testing such a simple 
loop. 

VI. 3. CONSOLE TESTING 

Testing with the aid of the control console requires some familiarity with the 
internal structure of the interpretive system (see Sec. XI) and with the basic language 
of the 650. Console testing is more attractive on the 650 than on most machines due 
to the ADDRESS STOP feature: If the CONTROL switch is turned to the ADDRESS 
STOP position, the execution of the program will proceed at electronic speed until 
the address set up on the ADDRESS SELECTION switches is reached. At that point, 
the machine stops, and the contents of various registers may be examined on the 
display lights or the program may be continued manually one step at a time. Alter- 
natively, the program may be punched out on cards at this point by merely feeding 
in one card with a PCH instruction, going into any vacant location, followed by a 
transfer card specifying this location. Console testing, in connection with the inter- 
pretive system, is likely to be needed only in exceptional cases. 

The ADDRESS STOP feature of the 650 may be used in conjunction with a spe- 
cial address stop transfer card when it is desired to start tracing from a certain 
instruction N in the middle of a program after running at full speed up to that point. 
(This may, of course, alternatively be accomplished using the tracing operations 
described in Sec. VI. 2, but then the value of N must be decided upon in advance 
and the proper program changes key punched. ) The procedure is as follows: Set 
the ADDRESS SELECTION switches to N and turn the CONTROL switch to ADDRESS 
STOP. Load as usual and run until the machine stops at the instruction N. (For 
details on possible earlier stops see below. ) Then set the CONTROL switch to 
RUN and load the tracing deck followed by the address stop transfer card, Tracing 
will begin immediately and the first instruction traced will be N. 

In choosing N it must be remembered that the loop box and location 000 are 
reset to zero when the tracing deck is loaded. If this restriction is inconvenient, 
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it can be circumvented by placing a special, card in front of the tracing deck. The 
card counter and the problem number are also reset to zero, unless the tracing 
deck has been modified to prevent it. 

If the CONTROL switch is kept in the ADDRESS STOP position when the pro- 
gram deck is loaded, one stop will occur when location N is reset by the memory 
reset card and another when the programmer's instruction is loaded into N. Also, 
stops may occur before instruction N is reached in the program, if N is referred 
to in an ADDRESS CHANGE or MOVE operation (but not if N is one address in a 
"conditional transfer instruction and control is transferred to the other address). 
After each stop, operation will resume when the PROGRAM START key is depressed. 
If the CONTROL switch is left in the ADDRESS STOP position during tracing, two 
stops will occur each time N is referred to (and one if N is the B-address of a 
transfer instruction). 

The program can be continued at full speed (punching mode) after a period of 
tracing by following the procedure described above with the punching deck in place 
of the tracing deck. 

The address stop transfer card has 69 1976 1952 24 1061 1098 in columns 
1-20 and a 12-punch in each of columns 1, 10 and 20. The special card for by- 
passing the resetting steps in loading has 69 1953 1952 24 1278 1953 70 1951 1344 
in columns 1-30 and a 12-punch in each of columns 1, 10, 20 and 30. (See Deck 7, 
Sec. XI. 7. ) 

If the value of N has been decided upon in time to get it key punched into a 
regular transfer card, (Sec. V. 3) this card may, of course, be used in place of 
the address stop transfer card in the procedure described above. 

A programmer familiar with the internal structure of the interpretive system 
will be able to think of many other cases where special needs can be met using 
maching language cards ("load cards 1 /). 
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VII SUMMARY OF OPERATION CODES 



l OPERATIONS 



2 OPERATIONS 







Page 






Page 






Page 


um. 


Alpha. 


Ref. 


Num. 


Alpha. 


Ref. 


Num. 


Alpha. 


Ref. 





GO to 


8 


000 


UNC STOP 


13 


300 


SQRT 










200 


COND STOP 13 


301 


EXP E 




1 


ADD 


10 


201 


TR SGN 


14 


302 


LOG E 




2 


SUB 


10 


202 


TR EXP 


14 


303 


SIN R ' 




3 


MPY 


10 


203 


TR 


14 


304 


COS R 




4 


DIV 


10 


204 


TR SUBR 


15 


305 


ART R 




5 


NGMPY 


10 


205 


TR OUT 


15 


350 
351 


ABS 
EXP 10 




6 


TR A 


19 


100 


LOOP A 


16 


352 


LOG 10 




7 


TRB 


20 


010 


LOOP B 


16 


353 


SIN D f 




8 


TR C 


20 


001 


LOOP C 


16 


354 


COS D 










110 


LOOP AB 


17 


355 


ART D 




9 


MOVE 


23 


101 


LOOP AC 


17 














011 


LOOP BC 


17 


400 


READ 


29 








111 


LOOP ABC 


17 


401 
410 


CONS 
PCH 


30 
26 








500 


SET A 


18 














050 


SET B 


18 


450 


START TR 


33 








005 


SET C 


18 


451 


STOP TR 


33 








600 


ADD A 


19 


452 


ST TR ERAS 


33 








060 


ADD B 


19 














006 


ADD C 


19 


454 


NOOP 


15 








700 


SUB A 


19 














070 


SUB B 


19 














007 


SUB C 


19 









36- 



VIII STOPS 



Error circuits in the 650 will stop the machine on attempts to use invalid 
information, such as, that represented by blank columns or double punches, as 
well as on several kinds of machine malfunctioning, and will indicate on the control 
console the nature of the error. If this occurs during the loading of a new deck, 
the cards should be examined. In other cases, a note should be made of the indica- 
tions on the console, and the procedure that led to the stop should, if possible, be 
repeated exactly in order to determine whether the error is systematic in nature. 

All stops, which are part of the interpretive system, will now be listed. 
Conditional stops will occur only if the PROGRAMMED STOP switch is set to 
STOP. On a conditional stop, the PROGRAM LIGHT in the OPERATING section 
of the console will be on and no lights in the CHECKING section should be on. The 
program will continue if the PROGRAM START key is depressed. On an uncon- 
ditional stop, the STORAGE SELECTION light in the CHECKING section will be on. 
Normally, operation should be discontinued after an unconditional stop and changes 
made in the program in order to avoid the stop. Alternatively, the program may 
be continued by having a transfer card (see LOADING) in the card reader, specifying 
the instruction to which control should proceed when the COMPUTER RESET and 
PROGRAM START keys are depressed. 

The location of the interpretive system instruction xxx on which the machine 
has stopped, may be determined by displaying the contents of location 1098 on the 
console. The display lights will show 60 Oxxx 1107. This process, called "monitor- 
ing u , may be performed as described in the 650 manual or, alternatively, by setting 
the storage entry switches to 60 Oxxx 8000 and depressing the COMPUTER RESET, 
PROGRAM START and PROGRAM STOP keys. 

If, in an exceptional case, it would be advisable to proceed to the next instruc- 
tion after an unexpected unconditional stop, this may be done manually as follows: 

(1) Set the CONTROL switch to MANUAL. 

(2) Check that the ADDRESS SELECTION switches are set to 1338. 

(3) Depress the COMPUTER RESET key. 

(4) Depress the TRANSFER key. 

(5) Set the CONTROL switch to RUN. 

(6) Depress the PROGRAM START key. 

As a result of this procedure, zero will be stored at C and 000 before the next instruc- 
tion is executed. If this is not desired, the ADDRESS SELECTION switches should 
be set to 1095 in step (2). To repeat the same instruction (on which the stop occurred) 
the switches are set to 1098. 

There is an alternative manual procedure for restarting after an unconditional 
stop which is simpler in the case of frequent use but is not recommended in general 
because it requires changing the setting of the STORAGE ENTRY switches. They 
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aroused in LOADING and must be set back to their normal positions for the next 

(1) nn t T t o h f fT S T T nnf AGE ENTRY SWITCH ES to 00 1951 1338+ (or 
vv 1VM 1095+ if zero is not to be stored or 00 1951 1098 + 
to repeat). 

(2) Depress the COMPUTER RESET key. 

(3) Depress the PROGRAM START key. 
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CONDITIONAL STOPS 



Address Lights 
1120 
1131 
1715 

1835 



Normal Cause 

Programmed COND STOP. (Display lights show B. ) 

CONS (Check STORAGE ENTRY switch setting. ) 

Loss of accuracy in SIN (Exponent of B exceeds 52) 
or COS. 

Loss of two or more digits of accuracy in LOG. 



UNCONDITIONAL STOPS 



2222 



r 

MOVE with C + A - 1 > 1000. 
•<PCH with B > C + 1. 



,READ with incorrect loc. or word count. 



3333 
4444 



DIV with B = 0. 
SQRT with B < 0. 



5555 



MPY with result out of range. 
DIV with result out of range. 



6666 



7777 



EXP with result out of range. 
LOG with B <_0. 

SIN with exp. of B exceeding 58. 
^COS with exp. of B exceeding 58, 



Cards missing or out of order in the system deck 
being loaded. 



9999 



Programmed UNC STOP (Display lights show B). 
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IX EXECUTION TIMES 

The execution times listed in this section are based on the standard 650 drum 
speed of 12, 500 r. p. m. They represent approximate theoretical estimates derived, 
in the case of the mathematical operations, from simple assumptions regarding the 
distribution of the numbers to be operated upon. For example, the part A ] of a 
floating -decimal number A = A l . 10 1 is assumed to be uniformly distributed 
between 1 and 10, although in physical problems there are reasons that favor a 
logarithmic distribution; extremely small and extremely large exponents are con- 
sidered very unlikely, etc. It is further assumed that the programmer has chosen 
storage locations on the drum without regard to timing, ignoring the fact that in 
the case of some operations the execution time will be a few milliseconds shorter 
for numbers stored in certain sections of memory. Some, but not nearly all, of 
the time estimates have been verified by tests. 

It should be stressed that the estimates of execution times are needed only 
for making comparisons or estimates of running time for problems or for choosing 
efficient ways of programming and will never affect the result of an operation. In 
comparing these estimates to estimates given for other interpretive systems or 
subroutines, it is important to verify by sample calculations or machine tests that 
the assumptions are realistic. 

To minimize the size of the table, the execution times listed refer to a basic 
case and corrections to be added in other cases are given at the beginning of the 
table. 5 
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650 INTERPRETIVE SYSTEM 
ESTIMATED AVERAGE EXECUTION TIMES IN MILLISECONDS 



(a) If A# 000, add 7. 2 ms. for ADD and SUB, 6. 3 ms. for MPY, NGMPY and 
DIV*. 

(b) If C ^ 000, add 6. 1 ms. for all mathematical operations, MOVE 000 and 
CONS*. 

(c) If the instruction has a minus sign, add 4. 8 ms. for all operations. 

(d) If, after a TR EXP or LOOP operation, control will proceed to the next 
instruction rather than to C, add 4. 8 ms. 



ADD 


65. 


7 


SUB 


65. 


7 


MPY 


67. 


2 


DIV 


74. 


3 


NGMPY 


67. 


2 



TR A 
TR B 
TR C 



37. 3 
37. 3 
42. 1 



MOVE 00037.7 
MOVE 40. 8 



(A = no. of 
words. ) 



UNC STOP 


28. 8 


COND STOP 


29. 8 


TR SGN 


19.2 


TR EXP 


24.0 


TR 


19.2 


TR SUBR 


44.4 


TR OUT 


26. 


LOOP A 


24. 


LOOP B 


28. 8 


LOOP C ( 


24. 


LOOP AB 


24. 


LOOP AC 


24. 


LOOP BC 


28.8 


LOOP ABC 


24. 


SET A 


55. 3 


SET B 


55. 3 


SET C 


55. 3 


ADD A 


44.9 


ADD B 


44.9 


ADD C 


44. 9 


SUB A 


44. 9 


SUB B 


44.9 


SUB C 


44.9 



SQRT 
EXP E 
LOG E 
SIN R 
COS R 
ART R 
ABS 
EXP 10 
LOG 10 
SIN D 
COS D 
ART D 

READ 



CONS 

PCH 



206 
197 
202 
192 
187 
238 

33.2 
187 
207 
240 
235 
271 

One card: 

101-4- 14n 

(n = no. of words. ) 

Succeeding cards: 

300 each. 

28. 8 

One card: 

163+ 12. 5n 

(n = no. of words. ) 

Succeeding cards: 

600 each. 



START TR 


28.8 


STOP TR 


24. 


ST TR ERAS 


38.9 


NOOP 


24.0 


TRACING 


600 per card 


LOADING 


300 per card 



*(See next page for footnote. ) 
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*Those who are particularly interested in time considerations may wish to know 
the exact increments on which the weighted averages in (a) and (b) are based: 

In ADD and SUB, 4.8 ms. if 17 < A < 41 (mod 50) 

9. 6 ms. if 1 <^ A< 16 or 42 < A < 49 (mod 50) 

In MPY, NGMPY and DIV, 4. 8 ms. if A> 17 or A =r 1 (mod 50) 

9.6 ms. if 2 < A < 16 (mod 50) 

In all mathematical operations, MOVE 000 and CONS: 

4. 8 ms. if 7 < C < 42 (mod 50) 

9. 6 ms. if 1 < C < 6 or 43 < C < 49 (mod 50) 

An easily remembered programming rule could be extracted from this in- 
formation: If locations between 17 and 41 (mod 50) are used for storing numbers, 
the increments given in (a) and (b) may be replaced by 4. 8 ms. 
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X SPECIAL TOPICS IN PROGRAMMING 



X. 1. SUBROUTINES; TRANSLATION 

A subroutine is a program expected to be of use as a part of the program in 
several problems or in several sections of the same problem. The mathematical 
operations in the interpretive system are indeed subroutines written in the basic 
language of the machine and reached through their operation codes, and anyone 
desirous of preparing an additional subroutine~of this type may avail himself of a 
vacant operation code (see Sec. XI) and write the program in a part of the memory 
below 1000. 

Subroutines written wholly or partly in the interpretive language may be 
reached conveniently using the TR SUBR operation defined in Sec. IV. 1. Suppose 
the subroutine begins at 900 and ends at 935. Instruction 900 should read: SET 
C 005 935 C 3and instruction 935 should read: TR 203 000C 1 Now suppose 
the subroutine is needed at step 700 in a program, and when it has been used, 
control is to be transferred to 680. Instruction 700 should read: TR SUBR 204 
680 900, The TR SUBR operation will take the quantity ("return address") 680, 
place it in the C-address of instruction 900 and then transfer control to 900. Instruc- 
tion 900, in turn, places 680 in the C-address of 935, and when instruction 935 is 
reached at the end of the subroutine, it transfers control to 680 as originally desired. 
Hence, the programmer using the subroutine only needs to know the identifying 
number 900; the transfer of control to and from the subroutine is handled by the 
TR SUBR in conjunction with the two instructions 900 and 935 provided by the writer 
of the subroutine. Subroutines needing only one input number and giving only one 
result (such as, the evaluation of one Bessel function for a given value of the argu- 
ment) will normally assume the input to be in 000 and will deliver the result there; 
in the case of several numbers, specified locations normally within the block 
occupied by the subroutine would be used for input and/or results. Subroutines may, 
of course, be used inside other subroutines without restriction. 

If the locations occupied by a subroutine are needed for another purpose, e. g. , 
another subroutine in the same problem, the subroutine may be translated to a 
different set of locations by a translating program developed by Miss D. C. Leagus. 
When the subroutine is written entirely in the interpretive system, the programmer 
is required only to separate data and constants from instructions, and the translating 
program will automatically decide which addresses of each instruction are subject 
to translation. Machine language instructions may also be used in a subroutine to 
be translated, provided certain conditions specified by the translating program are 
adhered to. 

Subroutines for the solution of cubic equations and of systems of linear 
equations have been written at the Laboratories. 
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X. 2. UNNORMALIZED INPUT; TRANSITION BETWEEN FLOATING- AND FIXED- 
DECIMAL FORM 

Nearly all of the mathematical operations in the system assume that the float- 
ing decimal numbers to be operated upon are in the normalized form defined in 
Sec. II. 2, i. e. , that the leading digit is different from zero unless the entire num- 
ber is zero. In processing empirical data, key punching is often facilitated by 
permitting leading zeros and reproducing a constant exponent. Such unnormalized 
data may be used in the interpretive system provided the first operation in which 
it is used is ADD or SUB with operand exponents differing by less than 10. 

A special case of unnormalized input is that of a zero with a non-zero machine 
exponent. If such a zero is added to a non-zero number with a smaller exponent, 
a number of digits equal to the difference between the exponents are lost. Con- 
sequently, zero should be equipped with exponent 00 unless the programmer knows 
in detail how the zero will be used in his program. Special provisions in the MPY 
and DIV routines make it possible to use a zero with machine exponent 00 in tnem 
without danger of exceeding the exponent range negatively. 

The converse problem of producing unnormalized output, e. g. , for the print- 
ing of tables in fixed-decimal form or for calculations in machine language is easily 
solved at the expense of one digit. Suppose for example that the numbers N. to be 
"unnormalized" or "fixed" are known to be less than 10 4 and output in the form 
XXXX.xxx is desired. Add the number 1000000054 (i. e. , 10, 000. ooo) to N if 
Nj > 0, subtract it from N. if N. < (using TR SGN) and punch. The output of the 
form + (10, 000 + |Ni|) , is ready to be printed on the tabulator with the leading 
1 and the constant exponent 54 suppressed by hammerlocks or wiring. If the num- 
bers are to be used in machine language, the 1 and 54 are shifted out. Rounding 
to a smaller number of digits is obtained by choosing the exponent of the additive 
constant (1000000054) correspondingly larger. 

X. 3. EXAMPLES 

In conclusion, two problems will be programmed in order to illustrate the use 
of many of the operations and methods described. 

First, suppose it is desired to evaluate the "error function", 

(1) *<»-£-£ e^t 

for a set of values x = a, a+& , a* 2 A, . . . , a+10A, using the RAND approximation 
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where 



(3) 



and 



ft~i/Cl+px) ; (p is a numerical constant)) 



(4) a^cx)^^ e~x* 

and to punch out the results as well as to store them for later use. The evaulation 
of the polynomial in Tl will be faster if (2) is written in the form. 



(5) 



^(x) = ^Cn(a l +TtCa^rL(a 3 4.n(a 4 4-n<a 5 )))))ce l (K> 



To make it possible to use the LOOP order in evaluating Gp*(x) this way, the 
coefficients a A will be stored in consecutive locations in decreasing order. The 

LOOP program will be given a form applicable to an arbitrary polynomial by 
including a "dummy" coefficient a = 0. Storage locations will be chosen as 
follows: 

Location Contents 

101-119 instructions (cards 1 -4) 

200 
201 

202 a I constants 

203 A f ( card 5 ) 
204 

221 x ""} temporary 

222 Q'U) r storage 

223 Tl J ("erasable") 

301 
302 

303 a 3 L coefficients 

304 a 2 \ in «£*C*) 

305 ai J (card 6) 
306 

401-410 Q*W J results 
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Card Loc, 



101 



The program might be written as follows: 
Alpha. Sign A or B 

SET C 4" 005 114 



102 



MOVE 



4- 



000 



1 103 


NGMPY 


4~ 


5 


000 


104 


EXP E 


4- 





301 


^105 


MPY 


4- 


3 


000 


^106 


MPY 


4- 


3 


204 


107 


ADD 


4- 


1 


000 


108 


DIV 


+ 


4 


201 


109 


MOVE 


4- 


9 


000 


110 


MPY 


+ 


3 


000 


rll-1 


ADD 


— 


1 


000 


| 112 


LOOP B 


+ 





010 


113 


NGMPY 


4- 


5 


000 


114 


ADD 


4- 


1 


000 


l 115 


ADD 


4- 


1 


221 


rii6 


ADD C 


4 





006 


' 117 


TR C 


4 


8 


114 


1 118 


PCH 


4 





410 


sJ19 


COND STOP 


4 





200 



202 

000 

000 

200 

221 

201 

000 

301 

223 

302 

005 

222 

201 

203 

114 

412 

401 

221 



500 



Next instruction in the problem. 



C Comments 

401 Set address of 
first <£*00 

221 First x is x = a 

000 -x 2 

000 e-x2 ^ J. 

222 (fc'OO'vSir*^ 

000 px 

000 1+ px 

223 H= f/C i+px) 

000 a 5 into 000 for LOOP 

000 prev. res.. 71 

000 add next coeff* 

110 loop in polynomial eval. 

000 polyn. 'Cf'00 

C 1 (J (x)= 1+ prev. res. 

221 x + A» next x 

001 next C£*(x) address 
103 test for end 

410 punch two cards 

500 end; stop, display last 
x, go to 500 on PRO- 
GRAM START 



An important remark should be made: If there is no shortage of storage 
locations and if the programmer does not mind writing a somewhat larger number 
of instructions, the running time for many problems can be decreased and the 
logic simplified by ,r unwinding M the innermost loop, i.e., by writing out the 
mathematical instructions in the loop in a straight sequence instead of using the 
LOOP operation. In the present problem, a sequence containing five MPY and four 
ADD instructions could replace the instructions 109-112 and also eliminate the use 
of the dummy coefficient a Q . The execution time for the polynomial loop would be 
reduced by nearly 1/3 and the LOOP operation could be used to replace the ad- 
dress change operations in the outer loop. The polynomial evaluation accounts for 
about 1/2 of the total running time of this problem. In many large problems, the 
innermost loop consumes an even larger fraction of the running time, making it 
important to program the innermost loop efficiently even at the expense of apparent 
inefficiencies elsewhere. 

The second illustrative problem reads as follows: For a given set of numbers 
x v , v =: 1, 2, . . . , 50, not necessarily equally spaced, the values of the Chebyshev 
polynomials T n (x v ), n=l, 2, . . . , 10, are to be computed using the recursion formula 
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(6) T T1+l Cx v )=ax v TnCx v )~T n -,Cx Y ), 

(T Q (x v ) ^ 1, Tj (x v )= x y ) and punched out in a compact form. 
In addition, the sum 

(7) f -ET.ofrv)] 2 (X v+| -Xy) , 

V = IV|-Xyl ' ' 

(X^ = 1) is to be punched out and the operator is to be given the option of also 
calling for the punching of partial sums of (7) at any time. 

Storage locations will be assigned as follows: 

Location Contents 



050 
051 
052 



0) 
1) 



constants (card 7) 



095-120 instructions (cards 1-6) 



199 
200 
201 
202 

210 

300 

301 

302 

400 4- v 
451 



The sum (7) and its partial sums' 
T Q (x v ) m 1 
Ti(x v ) 2 x v 
T 2 (x v ) 

Tl0<*v> 

2Xy 

[ T 10< x v>) 2 

[t 10 {x v I| 2 / v 



output 



1 



V 



temporary- 
storage 
("erasable") 



input 

(cards 8 - 18) 



In addition, locations 1-8 will be used in connection with a trick in program- 
ming the LOOP. 

The program may be written in many ways. The following is not necessarily 
the best: 
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A number of remarks are called for, many of them of general applicability: 

(a) The C-address of instruction 101 will, during execution, run through 
the values 000-008, but the result of the instruction is always called out from 000 
on step 102. This trick makes it possible to use the LOOP BC operation instead 
of address change, which is normally required if different sets of addresses are to 
be modified during a loop. 

(b) The instruction numbering was arrived at by starting the preparation 
of the program at instruction 101 with the intention of later adding an unknown 
number of initialization steps preceding it. This speaks in favor of not starting 
a program at 001. 

(c) The stop which would normally occur each time the CONS instruction 
is reached may be bypassed when found superfluous without any sacrifice by 
turning the PROGRAMMED STOP switch to RUN, since no COND STOP, SIN or 
LOG operations (the only other ones involving a conditional stop) are used. The 
operator decision regarding punching of partial sums is made using only the sign 
switch of the STORAGE ENTRY switches. This switch does not influence LOAD- 
ING. 

(d) The quantity x y is used so frequently that it was more economical to 
MOVE it into a fixed location than to apply address modification. The converse 
applies to x v +. j, which is used only once. 

(e) The constant 1 appears in three locations merely in order to simplify 
bookkeeping and loading, as well as changing the number of points x v in a later 
run. 

(f) An invaluable aid in determining whether the results of a calculation 

are correct is a mathematical identity which they must satisfy, and the programming 
of such checks is strongly recommended whenever it is possible. In the present 
problem, the identity 
1 

pr (xT) ^ 

(iii) \ 10 -* dx = X 



VT^ 



X 2 



2 



is closely connected with the computation of (7) if the x v are distributed over the 
interval (-1,1). 

(g) An alternative method of programming the outer loop, which would 
eliminate the address change operations at the expense of somewhat increased 
card preparation, would be to key punch the x v one to a card and give a READ 
order entering one x v at a time into a fixed location during the execution of the 
program. The difficulty arising from the need for x v + j on step 109 is not in- 
surmountable. 
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XI INTERNAL STRUCTURE OF THE SYSTEM 



XL 1. DETAILED DESIGN CONSIDERATIONS 

An expert examining the program at the end of this report will ask a number 
of questions about apparent duplication, about tight optimization in one routine in 
contrast to a lack of it in another, about the choice of operations and of methods 
of implementing them, etc. This section will attempt to answer some of these 
questions and also suggest a number of changes and additions that could be con- 
sidered for a second version of the system. Additional questions and suggestions 
from readers will be genuinely appreciated. 

In the early stages of system design, the following requirements were among 
those agreed upon, in addition to the general principles discussed in Sec. I. 2: 

(a) The arithmetic operations and those logical operations most likely to 
occur in inner loops (LOOP and certain TRANSFER operations) must be as fast 
as we know how to make them, regardless of the expense in storage. 

(b) The system must occupy at most 1000 memory locations. 

(c) The special functions must have full accuracy and unlimited range and 
most of them should be as fast as these requirements and available storage permit. 

(d) Optimum programming, (see the 650 Manual of June, 1955) in addition 
to being necessary for the attainment of (a) and (c), should be used locally in any 
program where the gain is significant but not at the expense of extensive rewriting 
of previously completed programs. 

(e) The programs must be so written that if the machine stops on any pro- 
gram step in a subroutine and control is transferred elsewhere before restarting, 
the subroutine, where the stop occurred, is left in a condition which assures correct 
operation the next time that subroutine is used. This implies that if a subroutine 

is used in more than one program, it must be initialized by each program rather 
than having a normal form used in one program and temporarily being changed at 
the beginning of other programs when needed there and then restored to normal at 
the end. 

(f) To facilitate changes, the individual programs (or ,r decks M , 1-20, see 
Sec. XL 7) that make up the system should be as independent of one another as they 
can be without excessive waste of^storage. This requirement was not fully adhered 
to near the end of the programming task. 

As a result of these requirements and of some oversights in programming, 
there are a number of storage registers which could be made available without any 
loss in system performance and a number which could be freed at some sacrifice. 
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A brief guide for finding such registers will now be given followed by a number of 
suggestions for their possible use in a revised version of the system. 

The 6 vacant C^-code locations and the 11 vacant registers listed in deck 5 
are, of course, available. The only distinction between them is one of mnemonics 
in connection with the choice of operation codes. In addition, it appears possible 
to salvage 22 registers essentially without loss by the following substitutions, but 
a careful check followed by machine testing is advisable: 



)eck 


Card 


Loc. 


Replace by 


20 


113 


1801 


1848 


16 


48 


1240 


1138 


18 


103 


1896 


1138 


6 


70 


1360 


1160 


16 


6 


1230 


1338 


18 


77 


1887 


1137 


2 


32 


1058 


1955 


12 


60 


1639 


1289 


16 


30 


1244 


1245 


19 


103 


1702 


1103 


17 


63 


1480 


1980 


10 


36 


1331 


1358 


8 


32 


1166 


1241 


12 


56 


1423 


1674 


17 


22 


1485 


1285 


18 


97 


1842 


1504 


17 


29 


1495 


1297 


5 


11 


1252 


1952 


5 


12 


1255 


1955 


5 


13 


1260 


I960 


5 


21 


1277 


1977 


5 


23 


1283 


1983 



It is, of course, necessary to determine, by sorting on instruction and data 
address, all places where the locations listed are referred to. 

Registers that may be freed at a price in speed include, above all, nearly 
40 extra registers used in the arithematic routines in calling out A and B, splitting 
them up and storing the parts. This is done separately in each of decks 12, 13 
and 14 to accommodate minor differences that facilitate optimization. To combine 
these steps without any loss of time is a task which, if possible, would require 
re -optimization of a substantial part of the system. At the expense of one revolution, 
they may be combined easily. Similarly, making the dissection of B common to 
all 02~routines would result in a substantial saving at the expense of lost time in 
cases (such as LOOP and TREXP) where B is irrelevant. To make this dissection 
common only to those routines where it is needed would be less profitable. 
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At some sacrifice in external characteristics, registers may, of course, be 
freed in any number of ways. If, in tracing, the modified instruction (redundant 
but convenient) is omitted, seven steps are eliminated. The MOVE operation for 
A# 000 is easily programmed in terms of LOOP BC and MOVE 000 and could be 
omitted, as could the special functions in degrees and to base 10 (or radians and 
base e, respectively). 

A number of suggestions for changes and additional operations will now be 
listed. Suggestions (1) - (3) use only the vacant registers and operation codes 
listed in deck 5 and can consequently be added to the system without difficulty at 
the option of any installation or individual programmer. For temporary use, they 
may be punched on separate cards and loaded after the system deck, in the case 
of (1) and (2) and after the punching deck, in the case of (3). Such cards should 
have an x-punch in column 5 and the four-digit location in columns 6-9. 

(1) Add an 2 -operation defined as follows: 

2 = 453 SWITCH Transfer on switch. Control is 

transferred to C if the Storage 
Entry Sign Switch is set to minus, 
to B if it is set to plus. 

This operation bypasses the stop that would occur if the same function were 
programmed by a CONS and a TR SGN order: It might be particularly useful in 
connection with tracing when it is desired to start tracing after a certain amount 
of running time has elapsed or for following the progress of a calculation by 
occasional punching of intermediate results at the discretion of the operator. 

The coding for SWITCH consists of the instruction: 



1453 10 



8000 1015 



Read console, 
routine. 



Go to TR SGN 



The execution time is 19. 2 ms. 

(2) Add an 2 -operation called COUNT having the same counting and testing 
properties as the LOOP orders but using a counter independent of the loop box and 
not capable of modifying instructions. Its function can be duplicated, e. g. , by a 
SET A, an ADD A and a TR A instruction. Its advantage lies in its speed and 
simplicity. The execution time is 24. ms. when control is transferred to C and 
33. 6 the last time when control proceeds ahead_and the counter is reset. A formal 
definition follows: 



2 = 800 



COUNT 



The number standing in the counter 
is increased by 1. Its new value 
is compared to B. If B is greater, 
control is transferred to C. Other- 
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wise, the counter is reset to 
zero and control proceeds to 
the next instruction. The 
counter is also reset in loading. 

If COUNT is used extensively, an expansion of the tracing program to punch 
out the contents of the counter, e. g. , in place of the problem number, would seem 
desirable. 



The coding for COUNT reads as follows: 



Loc. 



Op. 



1800 


10 


1314 


10 


1323 


11 


1381 


46 


1337 


10 


1396 


21 


1391 


16 


1066 


20 


1356 


00 


1317 


00 


1194 


20 


1378 


24 



Data 

1356 
1317 
8002 
1337 
8001 
1356 
8002 
1356 

0000 
1000 

0000 
1356 



Instr. 



Remarks 




0000 
0000 

1378*1 
1178J 



Call out and increase 

the contents, N, of the counter. 

Test N-f 1- B. 

On — , store N -+- 1 in the counter; go to C 

(in TR SGN program). 

On 4- , reset the counter, go to General 

Interpretation. 

The Counter 
Constant 

Change in LOADING to reset 
the counter. 



Note: If the COUNT program is loaded separately, the card loading 
zero into 1356 (step 9 in the program) must be included. 

(3) Include in the punching program (deck 9) a routine that prevents unwanted 
numbers from being punched out when the word count is less than the normal max- 
imum. This can be done on a tabulator with sufficient selector capacity (see Sec. 
XL 5). In the 650, it requires five locations and increases the execution time of 
the PCH order by 24. ms.when the word count is less than the normal maximum. 
The program, which may be punched on three cards, reads as follows: 



1949 


44 


1306 


1095 


1306 


20 


1980 


1307 


1307 


20 


1981 


1308 


1308 


20 


1982 


1309 


1309 


20 


1983 


1044 


1044 


20 


1984 


1973 



(4) Make room for the tracing program to be on the drum in parallel with 
the punch program, replacing the mode-of-operation deck (see LOADING) by an 
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ex- 



x-punch on the transfer card or a setting of the storage entry sign switch. An _. 
pansion of the loading program (about 7 steps) or of general interpretation would 
be needed, and the present overlap between tracing and punching is 34 registers. 

(5) If MOVE is omitted, except for A = 000, make this an 2 -operation and 
use the vacant Opcode for NGDIV. Alternatively, add a fast -operation, "NEG", 
identical with MOVE 000 except that it changes the sign of B. If NEG were available, 
however, it might be used in cases where, by slight reprogramming, a better pro- 
gram using NGMPY could be written. 

(6) Increase the number of logical operations, adding to the flexibility of the 
system and to the confusion of the beginner: Have a register called the "address 
counter", addressed, e.g., by 000 or by special operations and SET instructions 
referring to the address counter (as in 701 Speedcoding) where the present SET 
instructions refer to their own C-address. Have a set of TR A, TR B, TR C orders 
which automatically increase the address referred to by 1. These would have to be 
alternative to the address transfer orders in present use unless vacant 0, -codes are 
produced. * 

(7) Make use of addresses now ignored in some operations. For example, 
in CONS, use the B -address to call out a number B for console display when the 
machine stops. In TR OUT, make B a "return address" similar to that in TR SUBR. 
In START TR, or a new tracing order supplementing it, let B (or C) designate the 
number of steps to be traced before an automatic discontinuation of tracing. 

(8) Make Program Loading reset the registers below 1000 to zeros, unless 
told not to by an x-punch on the first card being loaded. 

(9) Have a conditional stop, or an operation effecting such a stop, on loss 

of accuracy in ADD and SUB, analogous to those in SIN and LOG. In many problems, 
particularly in connection with tests, such loss is legitimate, however,and a stop 
undesirable. 

(10) Replace or supplement the present error stops by the punching of an 
"error card". 

(11) Introduce an operation similar to ST TR ERAS, perhaps replacing it, 
which will cause the machine to trace the first, second and last repetition in a 
loop. 

(12) Add another LD-STD pair (at no loss in time) to General Interpretation 
(see cards 26 and 27) making ADD and SUB, as well as MPY and DIV available as 
internal subroutines. 

(13) Cut the execution time of several subroutines, such as, the arc tangent 
program, by making minor rearrangements, usually involving the expenditure of 

a few additional registers. 
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(14) Add an (^-operation, SPACE, which causes an x-punch to appear on the 
next card punched. 

(15) Interchange the functions of registers 1002 and 1702, causing the machine 
to stop sooner if a programmer accidentally attempts to continue upward from 
instruction 999. 

(16) Investigate whether a carry can ever occur on card 78, deck 18. If not, 
put the registers used on cards 78-82 to better use. 

(17) Replace or supplement the arc tan operation by an 0. -operation, ARG, 
which gives the argument (angle) of the point whose coordinates are (A, B). 

XI. 2. RELATED SYSTEMS 

Several systems supplementing the present one suggest themselves: (a) A 
system of symbolic or regional programming where the machine assigns absolute 
addresses in connection with loading; (b) A system externally identical with the 
present one, or very nearly so, operating on complex numbers, probably with real 
and imaginary parts in 8-2 floating -decimal form; (c) A system externally identical 
with the present one, or very nearly so, operating on double precision floating - 
decimal numbers, e. g. , 16-4; (d) A system of rr formula translation" or "automatic 
coding" (such as, the IBM Fortran for the 704) putting on the machine as much as 
possible of the burden of translation from a set of mathematical formulas to a 
program. 

XI. 3. NUMERICAL METHODS 

The study of numerical methods for calculation of the special functions included 
in the system was not nearly as exhaustive as would have been desirable and no 
claim to an optimal choice is made. 

The square root is computed by Newton's iteration method, 

where 1 <( B^ 10, using the initial approximation 
(2) x n slf .22B, . 



The evaluation of the trigonometric and exponential functions is based on RAND 
approximations (see Approximations for Digital Computers by Cecil Hastings, Jr. , 
Princeton University Press) to sin '& X and 10 x for < x < 1. Resembling the 
approximations obtainable by expansion in terms of orthogonal polynomials, these 
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approximations are in general somewhat more efficient than partial sums of Taylor 
series for a prescribed interval and accuracy, but it is not obvious that a further 
reduction of the argument followed by the use of a Taylor expansion could not have 
been better in the present case. For small x, as stated previously, the formula sin 
x as x (in radians) is used in order to retain significant figures. 

The logarithm and arc tangent are evaluated from fixed-length partial sums 
of power series after preliminary reductions of the argument, since eight-digit 
RAND approximations were not available. For the logarithm of B,, 1< B, <^ 10, 
the substitutions ~" 



O) U= jk ; V=± , for B,<e, 
(4> u = fe >V« 1.65, for B,£e, 



e 

(5) t =-a=i- . 
U+l > 

are followed by the evaluation of 

The constants, VlTand e 1 * ° 5 were arbitrarily chosen within the intervals that would 
lead to a minimal number of terms in (6). For x near 1, the logarithm is inherently 
less accurate than x since 

(7) d log X= ^ ^ iK 

and log x m whereas x^ 1. No substitution comparable to sin x = x can alleviate 
this difficulty. The use of a second (^-operation for log (1 -V- B) was considered 
but was rejected due to space limitations. This second logarithm could not be used 
to replace the present one for all values of the argument, since if the logarithm of 
a small number, say 10~ 10 is desired, the substitution, 

(8) I +- B = icf ,0 

forced upon the programmer, yields B= -1 exactly (in the eight-digit system used) 
with all digits of the input lost. 
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_ For the arc tangent, the reduction (after the argument is restricted to 0< 

B O by the obvious properties of the function) is based on the formula ~* 

(9) arc tan x = arc tan y +- arc tan %ZZ 

y l + Xy 

which is merely the addition theorem for the tangent rewritten. With y= .6, the 
use of (9) gives the desired accuracy in 

(10) arc tan z = Z fi - 1*+ f, *-,,.- Z?~\ 

L 3 S || J 

With Z ^ FXy f ° r x Z- 28 and z = x for x <.28. There is again some leeway in the 

choice of these constants. For small x, the substitution arc tan x = x is used to 
preserve significant digits. 

XI. 4. CONTROL PANEL WIRING FOR THE READ-PUNCH UNIT 

The control panel for the 533 Read-Punch Unit associated with the 650 is 
wired as follows: 

Col. 1, 1st Reading, to LOAD. 
R -h Sign, jackplugged. 
P-f- Sign, jackplugged. 

Col. 5, 1st Reading, to Pilot Sel. 1 X PU. 

Rd. Hold to PS1 Hold. 

Read Card C, Col. 6, to PS1 T. 

Read Impulse to PS1 N. 

PS1 C to Storage Entry C, Word 1, pos. 3 (from the left). 

Read Card C, Col. to Storage Entry C. 

7-9 Wd. 1, pos. 4-6 

10 Wd. 2, pos. 6 

H Wd. 3, Sign 

12-21 Wd. 3, pos. 1-10 

22 Wd. 4, Sign 



67-76 Wd. 8, pos. 1-10 

77-79 Wd. 9, pos. 4-6 
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Read Impulse 12 Wds. 1, 2, 9, 10, Sign 

Wd. 1, pos. 7-10 
Read Impulse ^ Wd. 2, pos. 7-10 

Wd. 9, pos. 7-10 



Word Size Emitter to Word Size Entry C 

10 Wds. 3-8, 10 

8 Wd. 1 

7 Wd. 9 

5 Wd. 2 



Storage Exit C to Punch Card C, Col, 

Wd. 10, pos. 3-6 1-4 

Wd. 1, pos. 3-6 6-9 

Wd. 2, pos. 6 10 

Wd. 3, Sign 11 

Wd. 3, pos. 1-10 12-21 

Wd. 4, Sign 22 



Wd. 8, pos. 1-10 67-76 

Wd. 9, pos. 4-6 77-79 

Wd. 2, pos. 10 ' 80 



Double Punch and Blank Column Detection as available and desired. 

XI. 5. CONTROL PANEL WIRING FOR THE TABULATOR 

The IBM accounting machine or tabulator used for printing from the cards 
associated with the interpretive system may be expected to perform some or all 
of the following tasks: 

(1) Automatic selection of different printing forms (i. e. , zero control and 
spacing between items on a line) for data cards and trace cards. 

(2) Selective spacing between lines. 

(3) Suppression of unwanted words from cards with word count less than 
5 (or 6). 

Since there are many tabulator models, it is not feasible to provide a detailed 
wiring diagram in this report. Instead, suggestions of general applicability will be 
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given. 

Exact selector requirements depend on the characteristics of each machine. 
As an example, requirements on a 416 will be given: The printing of signs requires 
6 single-position selectors with X-pickup and 6 positions of 11-12 separation (either 
special attachments or 6 positions of a selector transferred by an 11-1/2 impulse). 
Task (1) requires a digit selector (which may be put to duplicate use in task (3)) and 
34 selector positions with digit pickup (delayed pickup). Task (2) requires only 
3 one-position selectors on-the 416. Task (3) requires 55 selector positions with 
delayed pickup and some comparing units or a five-position selector and a digit 
emitter for control. 

The problem in connection with task (1) is to get the desired zero control 
and spacing in the two cases with the same setting of the hammersplit levers (also 
called zero suppression levers) and hammerlocks, on machines where zero control 
is not performed on the control panel. On a tabulator with 89 type bars, this may 
be done as follows: 



Type bar 



Direct Wiring 



Data Cards 



Trace Cards 



Alpha. 1-4 
5 
6-8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 



Col. 1-4, II*. 



Col. 7-9, II. 



Col. 5, II. 
Col. 6, II. 



Col. 11, I. 



Col. 12-21, II. 



Col. 22, I. 



Col. 23-29, II. 



Col. 11, I*. 
Col. 12, II. 
Emit 10. 
Col. 13, II. 
Col. 14, II. 
Col. 15, II. 
Emit 10. 
CoL 16, II. 
Col. 17, II. 
Col. 18, II. 
Emit 10. 
Col. 19, II. 
Col. 20, II. 
Col. 21, II. 
Emit 10. 
Col. 23, II. 
Emit 10. 
Col. 24, II. 
Col. 25, II. 
Col. 26, II. 
Emit 10. 
Col. 27, II. 
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Type bar 



Direct Wiring 



Data Cards 



Trace Cards 



Num. 



32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
1 
2 
3 
4 
5 
6 
7 
8-17 
18 
19 
20-29 
30 
31 
32-41 
42 
43-45 



Col. 30, II. 
Col. 31, II. 
Col. 32, II. 



Col. 41, II. 

Col. 42, II. 

Col. 43, II. 

Col. 44, I. 

Col. 45-54, II. 

Col. 55, I. 

Col. 56-65, II. 

Col. 66, I. 

Col. 67-76, II. 

Col. 77-79, II. 



Col. 33, I. 



Col. 34-40, II. 



Col. 28, II. 
Col. 29, II. 
Emit 10. 



Emit 10. 

Col. 35, II. 

Col. 36, II. 

Col. 37, II. 
Emit 10. 

Col. 38, II. 

Col. 39, II. 

Col. 40, II. 
Emit 10. 



(*) The symbol "11" denotes wiring from the second brushes (on some machines 
called "third reading") whereas u I ,r denotes wiring from the first brushes ("second 
reading") through a selector that separates ll's from 12 r s to the X-PU of an X- 
distributor ("pilot selector") through the transfer point of which an emitted 10 goes 
to the type bar in question. 

The hammer split (zero suppression) levers alpha. 4, 11, 25, 37 and num. 5, 
17, 29, 41 and the long hammerlocks alpha. 12, 16, 20, 24, 26, 30, 34, 38, 42 and 
num. 2 are raised. Left zero carry clips of width 3 are attached to hammersplit 
levers alpha. 6-8 and num. 43-45. 

Trace cards are distinguished by the presence of an 8 in column 80. This 
impulse is wired through a digit selector to the digit pickup of a row of selectors 
("class selectors" or "co-selectors with a controlling pilot selector") with a total 
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of 34 sets of points and to the hammerlock control hub. 

For performing task (2) in the manner specified in Sec-. V. 2. , four external 
switches are needed. On a 416, the single -double spacing lever and the minor, 
intermediate and major control switches can be used; on machines with four 
pluggable switches there is no problem. 

Spacing every 10 lines (alternative (b), Sec. V. 2. ) may be accomplished by 
adding 1 to a counter on card cycles and using the carry (which, on a 416, is 
automatically available at the counter total exit) to initiate a minor cycle during 
which spacing takes place and the counter is cleared. 

Spacing alternatives (c) and (d) both involve inspecting a card column at two 
reading stations and taking certain action when the digit at the first brushes is less 
than the digit at the second brushes. This may be done by wiring from second 
brushes to a comparing entry and from the corresponding comparing exit to the 
immediate pickup ( tr ZFS P. U. ,r on a 416) of a selector through the transfer point 
of which the digit at first brushes is passed. In case (d), this digit is taken directly 
to cause spacing; in case (c) it is wired to the digit pickup of a selector which 
initiates spacing on the next cycle. 

Task (3) is easily accomplished if a sufficient number of selectors are avail- 
able. Since the same task can be performed on the 650 at the expense of 5 locations, 
(see Sec. XL 1. ) the tabulator wiring will not be discussed here. 

XL 6. SYSTEM LOADING 

The interpretive system deck, normally with 6 words to a card and an x-punch 
in 5 to get the 1 in column 6 picked up as a leading digit of the address, is loaded by 
a deck of six self-loading cards (12 in col. 1) with 7 words to a card and a card 
number in the eighth word. The first card serves the sole purpose of making a 
fixed console setting possible. The System Loading program on these cards operates 
as follows: 



8000 
1951 
0004 
0053 
0106 
0153 
0204 
0056 
0055 
0155 
0102 



70 
70 
70 
70 
70 
70 
70 
60 
11 
44 
60 



1951 
0004 
0053 
0106 
0153 
0204 
0251 
0001 
0251 
0152 
0002 



133?\ 
0152 I 
0152 [ 
0152 f 
0152 ( 
OI52J 
0056 
005T) 
0155 Y 
0102J 
0007^ 



Read in the loading pro- 
gram from six load 
cards (B into 0001, C+ 1 
into 0002). 

Read a system card (non-load), 

Go to stop if expected loc. 
B 1 ^ loc. on card, L. 
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0007 
0205 
0057 
0051 
0052 
0103 
0156 
0206 
0054 
0203 
8002 
8003 
0151 
0201 
0207 
0101 
0105 
0107 
0104 
0157 
0005 
0003 
0152 

0001 
0002 
0006 
0202 
0154 



11 
11 
46 
21 
60 
10 
15 
10 
21 
11 
69 
24 
15 
10 
11 
44 
10 
60 
44 
60 
10 
21 
69 

00 
00 
24 
69 
00 



0001 
0252 
0152 
0080 
0006 
0251 
0202 
0252 
0070 
0252 
0)2531 

CLZJ 

0154 
8001 
0070 
0105 
0070 
0080 
0157 
0001 
0252 
0001 
7777 

B 
C 1 
0000 
0253 
0001 



0000 
0000 
0151 
8003 
0000 



Go to stop if C+ 1 <^ 
B' + n. 



Prepare accumulator for 
move. 

Store test constant. 



Move one word. 

Increase addresses by 1. 

Test for end of moving. 

Return to move another word. 
If C + 1 = B '-+- n, end of load- 
ing; go to console. 
Increase location by word 
count, go to read next card. 

Error stop. 



Constants 



A similar program is used for punching out the system in condensed form in 
case extensive changes, entered on self-loading, single-instruction cards, have 
been made. 

The Reset Memory Card, mentioned in Sec. V. 3. , is a load card with eight 
words. The program, essentially identical with one supplied by the IBM 650 Sales 
Research Group at Endicott, runs as follows: 



8000 


70 


1951 


1333 


1951 


69 


8000 


1953 


1953 


24 


0000 


1954 


1954 


69 


1957—" 


1955 


1955 


24 


0999 


1956 


1956 


61 


1958 


8003 


8003 


20 


[poof] 


0999 
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loo 



0999 

0000 70 1951 1333 



1952 8003 
0000 0000 



1952 00 0001 0000 

1957 11 1952 8003 

1958 20 0001 0999 

XL 7. PROGRAMS 

The complete programs of the system are listed on the next 21 pages, (i) - 

(xxi). In many cases, but not always, a constant used in two programs is listed 
in both. 
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650 INTERPRETIVE SYSTEM. 



DECK CARD LOC. 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 

1 
1 
1 

1 
1 
1 
1 



2 
2 
2 
2 

2 
2 
2 
2 

2 
2 
2 
2 



1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 



1 6 
17 

1 8 
1 9 
20 



2 6 
2 7 



1 
2 
3 
4 

5 
6 
7 
8 

9 

1 
2 



1095 
1014 
104 
8 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



06 
6 



80 



3 
1 
7 
1 
9 
6 
3 
1 
8 
5 
3 
7 
3 



1112 
1046 

1076 
1034 
8003 



21 1098 

22 1024 
2 3 1020 
24 1037 

2 5 1017 



17 29 
1735 



OR 

60 
10 
21 
60 
46 
30 
20 
60 
30 
20 
60 
30 
44 
10 
69 

1 1 
61 

60 
10 
65 

60 

69 
65 



69 
2 4 



1000 
1119 
80 02 
1096 

1200 
10 59 
8 002 
1016 



30 
1 5 
60 
69 



201 
015 
120 
069 



DATA 

109 
102 
109 
[ i+1 


102 
8 00 
000 
104 
8 00 
000 
106 
102 
8 002 

1017 
8 00 3 



INSTR. 



8 
4 
8 

] 

2 

3 
3 
3 
3 
5 
3 
3 
7 




101 
172 
800 
110 
106 
101 
02 
03 
04 
04 
10 
06 
07 
800 
108 



104 6 
106 1 



1. GENERAL INTERPRETATION . 

Increase i to i+1. 
(see cards 26 & 27 below) 
Start dissection system instr. 
Shall the loop box be used? 

Store C* 



Store B* 

Is an 02-operation called for? 

Go to an 0| -subroutine „ 

Add the loop box to the instr. 
Remove the minus sign. 



8001 1034 
1037 8003 
10 4 5 [IOOO+O2] Go to an 02~subroutine, 



I A "i 

0001 

8 002 
104 5 



1732 
153 9 



28 17 32 35 0002 



110 7] 

V Constants 
108 1 

1 00 OJ 

Loop box (initially 0) 

17 3 5^ Restore the multiplication 
^■9A^ \ routine (Oj-3) to normal 
(The trig, program makes 
1 4 4 5 J special use of it). 



W L 



00 3 
10 30 

9 99 9 

000 3 
116 2 
[ B ] 

10 9 
1120 
102 3 
8 2 



11191 
800 2 
109 6. 
112 

10 5 9"] 
800 2 
101 6J 
112 



2 # TRANSFER OPERATIONS . 

OyOOQ . UNC STOP 

Call out Hi "for console display 

Stop. 

Q£-20o7 q3Nr) stop 

Call out B for console display 
Stop; Go to C on PROGRAM START 



!J 



0a°201 . TR SGN 



10 15 __ 

1 6 9> Test the sign"of prer. result 

10 2 7 Go to C If -, 
10 2 7 Go to B if +. 



2 



13 1202 60 1009 1013 



02-202 . TR EXP 



-ii~ 



2 
2 
2 
2 
2 

2 
2 
2 

2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 

2 
2 
2 

2 

2 
2 
2 
2 
2 



3 
3 
3 
3 
3 



19 
20 
21 



2 2 
23 
24 
2 5 
26 
27 
2 8 

2 9 

3 
31 
32 

33 
34 
3 5 



37 
3 8 
39 
40 
41 



1 
2 
3 
4 
5 



1013 
1021 
1029 
1036 
1012 

12 3 
1027 
1035 



20 
12 
04 
15 
16 
800 
107 
13 
09 
04 



058 



120 5 
112 8 
8003 



3 6 1454 



1 001 
1064 
1071 
10 7 7 
1135 



30 
67 
30 
16 
4 6 

6 
30 
10 



65 
35 



67 
21 
60 
60 
60 



30 
69 
30 
10 
10 



0002 
8002 
0001 
104 5 
109 5 

102 3 
000 3 

103 8 



102 
000 
114 
105 
111 
[ c 
000 
8 00 
104 
000 

[ c 



102 3 
000 3 



00 5 
1017 
0002 
8 001 
113 8 



102 1 
102 9 
10 3 6. 
1012 
1120 



Get the exp. of prev. result 

Compare it to B; Go ahead 

if exp. < B, to C if exp. ^> B, 



10 2 7") 02-203 . TR 

10 3 5V Replace i+1 by C in 

10 4 3 J general interpretation. 



112 71 
104 
115 1 
61 



02-204 . TR SUBR 



Set the C-address of the 
Y instruction at C equal to B; 
So to TR to C. 



11 

800 2 
1075 
113 3 
109 2 
1049 
105 8. 
12 3 J 



8 3 V Go to C in machine language 
[ C ] J 



10 9 5 




10 6 
107 
107 
113 
109 



Og-454 . N00P 

Go to general interpretation < 



Constants 



3« LOOP OPERATIONS 
02*001 . LOOP C 
Move B to the C-address 
" position; Call out and 
increase the loop box. 



3 
3 
3 
3 
3 
3 

3 
3 
3 

3 
3 
3 

3 
3 



6 
7 
8 
9 
10 
11 



093 
051 

54 
062 

1 55 
113 



1 8 
19 



1100 
10 57 
10 2 5 

1008 
1074 

1110 
1121 



3 
10 
10 

30 
10 
10 

69 
30 



8 002 
105 4 
8001 
1017 
8 002 
1017 

0001 
1017 
10 3 1 

00 1 
1017 
1031 

1017 
0001 



10 5 1] Common LOOP steps 

115 5)" Compare the loop box to B; 

10 6 2 if less, loop again (go to C). 
112 0J 

1 1 1 31 If equal, reset zhe loop box, 
1 9 5 J go to general interpretation 



10 5 7 
102 5 
1 09 3 

100 8 
1074 
113 5 

112 1 
102 8 



Oa-100. 


LOOP A 


02-101, 


LOOP AC 


Og. 


•110, 


LOOP AB 
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3 
3 

3 
3 
3 
3 

3 
3 
3 

3 
3 

3 
3 
3 
3 
3 



4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 
4 
4 
4 
4 

4 
4 
4 
4 
4 
4 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 



20 
21 



1028 
1136 



22 1111 

2 3 117 1 

2 4 1078 

2 5 118 5 

26 1010 

27 13 43 

28 1003 



29 
30 

31 
3 2 
33 
3 4 
3 5 



12 
13 
1 4 
1 5 
16 



7 
8 
9 

1 
9 



23 
2 4 
2 5 
2 6 
27 
2 8 

2 9 
30 
31 

3 2 
3 3 



1011 

118 



138 
031 
47 
039 
139 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 



09 



153 
312 
322 
3 30 

3 88 



67 
28 
40 
63 
72 
36 



1197 
1176 
80 03 
8002 



10 
10 

69 
30 
10 
10 

30 
10 
10 

30 
10 



65 
4 5 
30 
16 
20 
3 5 



10 
4 6 

10 
30 
10 
16 
1 5 




6 
1 

1 
6 



8001 

103 9 

1017 
0001 
8001 
113 9 

000 4 

104 7 
1017 

000 4 
113 8 



0100 



01 
01 



00 
00 



3 
4 
3 
3 
3 
5 



8Q02 
110 9 
000 
10 2 
128 
000 
102 
131 
8 001 
104 5 
136 7 

8002 
000 3 
132 5 
128 3 
114 2 



102 
8 00 
104 
00 
132 
114 




113 6 
109 3 

117 1 
107 8 

118 5 
109 3 

13 4 3 
1003 
109 3 

118 
13 4 3 

000 1^ 

100 
1000 
100 1J 



02*111 . LOOP ABC 



02"010 . LOOP B 



02-011 . LOOP BC 



Constants, 



100 
118 
171 
132 
1.18 
109 
132 
113 
119 
109 
115 



\ 



4. MOVE 
Ql-9 . MOVE 
Test N (i.e 



Store N-l 



1312' 

132 2 

133 
138 8 
119 7 

132 8* 

134 
136 3 
117 2 

133 6 
119 7, 

1176 
800 3 
80 2) 

1102} 




• 9 A) | go to special 
move if N-0 



Go to stop if N-l+C £ lOOOo 

Test for upward or downward move. 

Initialize for downward move. 



Initialize for upward move. 



Store test constant * (see next p.) 
Complete initialization 

Move one word 

Test for completion 

Change addresses by 1, 
return to move next word 
If N-0, Get E and go 
to storing routine. 
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4 
4 
4 

4 
4 
4 



5 
5 
5 
5 
5 
5 

5 
5 
5 

5 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 

5 
$ 
5 
5 
5 
5 
5 
5 
5 
5 
5 



6 
6 



1 

2 
3 
4 
5 
6 

7 

8 
9 


1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 


1 
2 
3 
4 
5 
6 
7 
8 
9 




1363 
80 02 
1315 

1325 
114 2 
1065 



06 
07 
08 
09 
53 
00 



1009 
1023 
1 4 5 



1 

1 

1 

1 

1 

1 

1 

1 

1-2 

1 2 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



5 
5 
5 
6 
6 
6 
6 
268 
70 
7 
7 
7 
7 
8 
8 



2 
4 
7 
8 
3 
5 



89 
91 
93 
94 



297 



35 
31 
31 
32 
37 
38 
33 
39 
39 
6 



10 4 4 

10 87 

10 8 8 



15 106 5 
60 [ B ] 
69 2222 

69 
24 
60 



77 
77 
77 
77 
77 



88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
8 8 
88 
88 
88 



69 
69 



77 
77 
7 7 
77 
77 



55 5 55 5 
55 5 55 5 
55 5 55 5 



8 88 8 

8 88 8 

8 88 8 

8 88 8 

888 8 

8 88 8 

8 88 8 

8 8 8 8 

8 8 8 8 

8 88 8 

8 88 8 

8 8 8 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 



7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 



7 7 
77 
77 
77 
77 
77 
77 
77 
7 7 
77 
77 



80 2( 
14 4 5) 
10 9 5 Error stop for MOVE, READ & PCH. 



:} 



80 2" 

110 2V Constants 
144 5. 



77771 
777 7 
7777 
777 7 
777 7 



5. ERASABLE AND VACANT LOCATIONS , 



r Vacant 02-codes 



77 7777 777 7J 



1 1 4 
1341 



55 5 5 
555 5 
555 5 

88 8 8\ 

88 8 8 

88£ 8 

8 8 8 8 

888 8 

88 8 8 

888 8 

888 8 

888 8 

888 8 

88 8 8 

888 8 

88 8 8 

8 88 8 

8 8 8 8 

888 8 

8888 

888 8 

8 8 8 8 > Erasable 

88 8 8J 



Inter-subroutine storage 



INITIALIZATION CONSTANTS 
FOR MOVE 




Downward 


Upward 


# 


B+N-l 


-B 


## 


B 


-(B+N-l) 


■JBtt* 


C 


-((MM.) 



777 

777 

777 

77 7 

77 

77 

77 

77 

77 

77 

7 7 



7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7J 



> Vacant registers 



6. ADDRESS CHANGE OPERATIONS. 



119 5 



119 5 0^, TE B 



_v- 



6 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 
6 
6 
6 

6 



3 1089 69 1342 1195 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 
6 
6 
6 
6 



4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 

6 
7 
8 
9 



2 



21 



2 3 

2 4 

2 5 

2 6 
2' 7 

2 8 

2 9 



3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
4 
4 
4 
4 
4 

4 
4 
4 

4 
4 




1 
2 
3 
4 
5 
6 
7 
8 
9 

1 

a 

3 
4 

5 
6 
7 
8 
9 



1195 
10 53 
1311 
1369 
8 02 
1198 



59 
17 

25 
83 
92 
99 



1140 
1341 
1342 
1072 

100 5 

10 50 



63 
69 
77 
8 4 
292 
349 
1 57 
164 
170 



8 02 
12 89 
1108 
1122 
10 32 
1144 



1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 



1 58 
115 
118 
1 30 
196 
274 
267 



103 
116 
8 
16 
104 



5 1006 
5 1 10 07 



65 



2 2 1500 65 



20 
60 



3 5 
3 

21 



65 
66 



8 
2 
3 
2 
1 



119 
8 00 
000 
107 

[ A 

[1,4,7] 
8 002 
000 3 
8 00 
00 
104 
120 



00 01 
00 4 
00 7 



6 5 8 3 



1103 
110 4 

1116 

127 4 
108 

128 9 
104 5 

00 3 
1160 
126 7 
107 3 

1 B 3 
[0,3,6J 

126 5 

8 3 

00 3 

8 00 3 

1115 

10 2 3 

10 2 3 

3 

126 5 

[0,3,6] 

[ B ] 

3 



00 6 



Common steps . 

Set amount of shift for TR A, B or C. 



Get the 
located 



instruction 
at A. 



go 



Separate out its A-, B-, 
or C-address. 



ii 



Compare this to B. 
Go on if equal, to C 



Constants 



if unequal. 



116 3 Og-005 . SET C 
116 3 02*050 . SET B 
116 3 0p«500 . SET A 



9*1 



) 



Comm on steps 

Set amounts of shift. 



Set the address for storing 
the modified instr. at B. 
Get the instruction (!) 
located at B. 

Store the right end of !. 

Destroy the old A~, B~ or 
C-address 



09 



' Replace it by C 



Attach the right end; 
store at B. 



12 67 

110 8V Constants 

109 5 



1023 1079 02O06 . ADD C 
1023 1079 Qp-007 . SUB C 



-vi- 



5 2 1079 30 0007 1004 Common to both 



6 5 3 1060 

6 54 1070 

6 5 5 1129 

6 5 6 1600 

6 5 7 1700 

6 58 1189 



65 1023 1129 02*060 . ADD B 

6 6 1023 1129 02»070 . SUB B 

30 0004 1004 Common to both 

65 1023 1189 Og-600 . ADD A 



66 102 3 
30 0001 



118 9 02"700 . SUB A 
10 4 Common to both 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 

^6 
6 



70 
7 1 
7 2 



00 
6 
39 
3 5 
31 
37 
8 00 
134 
110 
115 
126 



1360 
1123 
1073 



21 
60 
60 



126 
104 5 
3 
1360 
26 7 
12 3 

B 

15 

26 
26 

B ] 



1 
1 



1156 



[ 








jap. 



12 6 
109 



109 5' 
134 8 
12 8 9 



Common to all Agp & SUB 
Store shifted ±C; " 
Set address for storing 
modified instr* at 8 

Get 1 



Add ±C to j Bf | store 
(with original sign) at B # 



Constants 



7 
7 
7 
7 
7 
7 
7 
7 

7 
7 
7 
7 
7 
7 
7 
7 



1 
2 
3 
4 
5 
6 
7 
8 

9 

1 
2 

3 
4 



8 00 
1333 



78 
75 
44 
4 2 

77 

76 



70 
60 
24 
20 
24 
2 4 
60 
4 5 



15 1038 60 



951 
242 
27 8 

017 
98 6 
952 
136 4 



95 9 
98 5 
97 6 
06 1 
951 
38 



8 
8 


1 

2 


1401 
1131 


01 
6 


8 000 


8 
8 


3 
4 


14 00 
1174 


69 

24 


1318 

1 2 i 8 



133 
112 
119 
117 
117 
134 
104 
137 



} 






7. PROGRAM LOADING , 

Read the first program card. 

Set an instruction for return 

from the READ routine 

Reset 000 

Reset the Loop box 

Reset the card counter 

Test the word count, n. If V 0, 

go to READ to move n words 



17 1242 70 1951 1344 



1 3 7 6\ If n*0, put the problem no, 

1 3 9 7 J into the punch band, 

13 4 5") decide if tracing is required, 

1 1 8 J (dep, on contents of 1976*), 

112 61 store the address of the 

17 2 9 > first programmer instruction, 

go to execute it (gen. int.) 
1 1 7 J f30 0003 1019 if not tracing 
21 1980 1386 if tracing 

Constant (return from READ) 



#1. &(> f>.y 



« 



8. READING OPERATIONS 



1 *• 3 1] 0g»401 . READ CONS 

144 5) Stop* On start „ read console 

switches, go to storing routine, 
117 4] 02=400 . READ 
1 4 8 6r Set return instructions, go to 



•vii- 



8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 
8 
8 
8 



5 1486 6 9 124 1 116 8J steps common with PC H 



9 
9 
9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



go 

2 1 



2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
2 



2 
3 
4 
5 
6 
7 

a 

9 

o 
l 

2 
3 
4 
5 
6 
7 
8 
9 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



272 
736 
326 
365 
7 33 
347 
372 
382 
740 
590 
3 64 
6 86 
56 
231 
09 4 
8 02 
8 03 
13 20 
335 
193 
37 3 
339 
2 78 
39 
346 
584 
166 




3 3 1188 

34 1318 

3 5 1241 

3 6 1324 



69 
60 
21 
2 4 



195 3 
128 5 
127 7 



80 3' 
139 
173 6 
132 0. 



Return from common steps* 
Read a card. 

Go to stop if loc. L on the card 
differs from progr. first loc« 
(the stop is in MOVE) 
Go to stop if B f +n > Ol 
(n=word count, B f "current 
first location) 



Calculate and store 
test constant 

Complete initialization 

Move one word from the read 
band into place 
Increase addresses by 1» 

Test for end of moving« 

Return to move next word 

If B'+n * C+l, end of READ 

Add n to B* to get new B', 
go to read next card* 

Constants 



1 1410 69 1964 1168 



9« PUNCHING , 
02»/qO a PCH 



1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



168 
81 
2 32 
361 
319 
134 



2 4 
65 
10 
30 
1 5 
20 



72 
89 
98 
66 
73 
88 
72 
62 
99 
61 
68 
48 
7 



127 2 
10 2 3 
104 5 
3 
10 2 4 
1291 

197 7 
8 2 
196 6 
196 9 



10 8 1") Common with READ 

12 3 2 * ~~ 

1 3 6 1 1 Prepare for testing 
1 3 1 9 f and initialization. 
113 4 
127 2) 



9 

2 

9 

9 

97 
8 
00 
196 



197 4 



98 
9 9 
31 



8 



1 



3 

8 

21 

2 

9 

1) 

81 

8 



0, 



PCH program only 

Go to stop (in MOVE) if 

B > C+l. 

Set the word count 

Let B r -first loc, not yet 

punched. If B'+5 < C+l, 

go to punch 5 to a~"card. 

If not, go to further testing. 

Prepare to move n words 
to the punch/band 



-viii- 

9 21 1950 15 1977 1987 

o o5 H 87 10 1990 8002 

o is «nn§ § 9 I B ' 1 8003> Move one word 

9 34 8003 2 4 [1 9 7 91 [1990+nl] 

9 *§?£ 6 L 19 8 6 J L 19 7 1 ,J Gotopunch 

9 2 6 1993 15 1024 198 71 

o IZ H 93 15 102 <* 1987/ 

o §g 1994 15 1024 

9 2 9 1995 15 1024 

9 30 1996 15 1024 

9 31 197 1 10 l_Q-2 4 

9 32 1967 21 1986 

9 3 3 1997 71 1977 

9 34 1947 60 1977 



Return to move 
another word 

Increase the card no. 

Punch a card 



9 If if«f I? HZS 19 6 5 Uet the location no. 
9 3 7 ? £ I ? 11 ! 9 6 SJ for the next card 

9 38 1949 44 f 9 2 3 1 ^ 4 9| If B -C + l, punching 
jo A * 4 y 44 1973 109 5Jia completed. 

9 39 197fi 30 n n n •* iaioi J?* 8 const « kee P s tra e« orders 
9 40 19 64 2? 1077 {inoperative during non-tracing. 

9 41 19 69 oooZ ! 9 8 9 (This const. «ill be 00 0000 1095 

9 42 199 99 9 9tes^r Cin8(FCHta0peratiVe) 

y 43 1963 24 1978 199 1/ 

9 45 1958 12 1Q7P ?S 3 (When the PCH progr. is loaded 

40 XJbB 34 1978 199 lj] ^th jEir^gram loading, this instr. 

Lwill load into 1963 (!). 

•in *i -irv^-i /-, ., .. 10 ♦ TRACIM G 

10 2 ?n nnn^ 1386 ^^iioT^f ^en. int. 

10 % ^Ist on ? ° ° 3 1946 Store ^odifiedinstrT for trac. 

Iq 3 * ? 4 £ 20 1023 1947 Store C 

10 
10 
10 
10 
10 
10 

10 ta 19^ SS 22SS 1966^CalloutIifO,/0 
10 

10 

10 16 J52S ^ 1Z°A J 9 ^ Store A or Of or tracing 

10 

10 1§ 80o£ ^ ?l", BO 2^StoreI for tracing 
10 

10 32 1974 34 1984 ^ ^ } Store P rev - re ^t for tracing 

1U 241987241981198 8J 



10 
10 




1 9 8 91 Store instr. no. for tracing 

199 0. ~ I j)L , v , i.„ 



,yww •■■■.- 10 & r x^yu 15 1997 8002 

^ " 10 I 9 1 9 9 1 ■ 1 ? £ ^ 9 f § M5^ °^ inal *-f • *" tracing 

10 30 1992 65 19 8 6 199 31 



*ix- 



k&(f 



1 
10 
10 
10 

10 
10 

10 
10 

10 
10 
10 
10 
10 
10 

10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 



11 



12 
12 
12 
12 
12 



3 1 1993 

3 2 1994 

3 3 199 5 

3 4 1996 



3 5 
3 6 

37 
3 8 



4 5 

4 6 

4 7 

4 8 

4 9 



5 
5 
5 
5 
5 
5 



2 
3 
4 
5 
6 
7 



1 
2 

3 
4 



6 
7 
8 
9 
10 



14 50 
1331 

14 51 
13 58 

1452 
1332 
1091 
1165 
1334 
8 02 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



99 7 
998 
999 
141 
173 
1 54 
9 78 
9 64 
976 
24 
9 58 



4 04 

414 
42 7 
4 30 
4 39 
441 



15 
2 
7 1 
60 

69 

2 4 

6 9 
2 4 

69 
24 
6 5 
16 
69 
2 4 




6 
6 


3 
3 



21 



4 4 



1 14 45 21 



21 
65 
4 5 
30 

1 5 
6 9 



102 4 
198 6 
19 7 7 
1260 

197 6 
1061 

115 4 
1061 

197 6 
106 1 
109 8 
1173 
1141 
[i+1 } 



5 400 

00 3 
000 6 

198 
1, 
196 5 



100 
102 
14 3 

144 



880 2 24 [C] 

9 1443 2 4 



10 82 
149 2 
14 49 
14 57 



24 
65 
30 
1 5 



128 9 
104 5 
3 
146 



199 4> 
199 5) 
199 6 
110 5 



Store card no. for trace cards 

Punch a trace card 
Return to normal gen, int. 



1 3 3 1) 02<50 . START TR 

1095)" Modify gen. int. to include 



tracing expansion 
13 5 81 Qg-451 . STOP TR 
109 5/ Restore gen. int. to normal 



1332] 02-452 . ST TR ERAS 

10 9 1} Modify gen. int. to include 

1 1 6 51 tracing expansion. 

1334) 

8 2V Replace the present progr. 

1 9 5 J instr. by a NOOP 



08 8 
196 
197 

001 
101 
000 
109 
1 38 



4) 

8 

2 



196 8 J 




r Constants 



Word count }£ trace identif . 
Make PCH inoperative when tracing. 
Will be 30 0003 1019 when not trac. 
When the trace program is 
loaded with program loading, 
this will load into 1963 (l). 
11. STORING THE RESULT 
[Common to all math, routines ) 



Store C in 000 for prev. res. and 
in 1009 for condo transfer ops. 
If C-000, go directly to gen. int. 

If C / 000, store C iii C, 
go to gen. int. 



10 9 5 Constant 



149 2 

144 9 

14 5 7 

800 2 



12. ADDITION AND SUBTRACTION 
0j»l , ADD 

Get B 



5 8002 60 [ B] 147 3. 



10 8 3 
15 4 2 
1599 
120 7 
8 02 



11 14 7 3 

12 148 3 

13 14 93 



24 
65 
30 
1 5 
61 

21 
6 5 
4 5 



128 9 
10 4 5 
3 
1310 
[ B ] 

127 8 

128 9 
149 6 



15 4 2 0i*2 . SUB 

159 9) "^ 

1 2 7 I Get -B 
800 2 
147 3. 



148 3 

149 3 
149 7 



Common steps 

If A-OQQ, get TL directly, 



-X- 



1 
1 
1 
1 
1 
1 

2 
2 1 
2 2 
2 3 
-9-4 
2 5 
2 6 
2 7 
2 8 
2 9 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 




1 
2 
3 
4 
5 
6 
7 
8 
9 



40 
41 



149 
14 
8 
14 4 

4 

4 

4 

4 

4 

5 

5 

6 

6 

6 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



9 

1 
6 

7 

8 

8 

9 



6 

1 
63 
43 
4 4 
24 
408 
429 
43 
48 
49 
65 
46 
41 

7 



1 

4 

9 



2 

7 



48 
5 4 



42 
4 6 



00 3 
140 9 
[ A ] 



2 
21 



7 

8 

9 
2 
5 
5 
5 
5 
4 
1 
6 
6 
8 
5 

5 
7 
5 

6 
1 


6 
9 

5 
6 


6 
5 
9 
1 
5 
9 
4 
5 
9 



00 
126 

1 27 
1 27 
000 
1 29 
125 
8 00 
126 
4 
143 3 

63 

69 

25 

27 
00 
148 

165 
145 
126 

8 
145 
126 
149 
8 00 



8 02 
8 3 
0002 
1461 



000 
141 
142 
8 00 
127 
147 
125 



127 2 



c 

8 00 
129 
126 
127 
125 



140 
800 
144 

4 

4 

4 

4 

4 

5 

5 

6 



60 

66 



40 
5 2 
4 7 
48 
54 



2 V Get A if A ^ 000 
6J 

5 
5 

11 - 

8 V A - Ai ,a # Split up 

5 J and store 
3] _ 

9 I B - Bi,b, Split up 

6 J and store 
3J 

OYPut |b-a| into shift instr. 

2J 
71 

8 I If a*b, form Ai+Bt«CJ„ 

9 

5J 

°>I* l c i| > 10 » c i " c '/io. 

3 ] 

If If Ci < 0, get c«a+l, 

1 r combine with Cf , 

7 go to store. 
5J 

ll Same if C^ > 

1) 

5 If Ci»0, go to store machine 0. 

y 
1 

9 
7 
4 
1J 

3 If a / b, prepare to shift* 
5 

3 

5 

5J 

3 

21 

9 

3 

9 

1 

9 1 If b > a, interchange 

2 f Ai and B-\, replace a by b. 
9J 



r If |C{ | < 10, get Ct by 
shifting! correct exp. 



„ If_|b-aj > 10, get OA or 
OB, go to store. 



- If b < a, shift Bt • 



JL (i 

12 
12 
12 



6 6 
6 7 
6 8 
6 9 
70 



1310 
14 6 
14 9 
1433 
1456 



61 
6 
60 
31 
1 



147 3" 

147 3 

X 4 4 6y Constants 

1413 



13 
13 
13 
13 
13 

13 
13 
1 3 
13 
13 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 



1 
2 
3 
4 
5 

6 
7 
8 
9 
1 



1 

1 
1 

1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 
2 
2 



2 8 
29 



3 
3 
3 
3 
3 
3 

3 
3 
3 
3 
4 
4 
4 
4 
4 
*4 
4 
4 
4 
4 
5 
5 
5 
5 




1 

2 
3 
4 
5 

6 

7 

8 

9 



1. 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 



10 84 
1592 
15 49 
15 57 
8 02 

10 86 
1692 
1649 
1657 
8 02 



5 78 

5 86 

59 

59 

60 
8 
159 



1 

1 
1 
1 
1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
-1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



88 
9 1 
4 6 
64 
4 8 
59 



717 
27 3 
6 77 

5 33 

6 28 

4 84 

7 43 
9 9 
54 
7 
58 
53 
63 

571 

5 28 
581 
5 32 
539 



24 
65 
30 
1 5 
60 

24 
65 
30 
1 5 
61 



20 
60 
30 
4 4 
2 1 
68 
15 
1 8 
20 
4 6 
4 4 

9 
5 
5 



3 5 

4 4 
69 



128 9 
104 5 
00 3 
1560 
[ B ] 

128 9 
104 5 
00 3 
166 

[ B ] 

128 3 

128 9 

159 6 

00 3 

150 8 

[ A J 



2 
5 

7 
3 
2 
7. 
2 
2 

56 8 

27 7 

28 

59 

59 

2 5 

27 
00 4 
8 3 



00 

1 51 
127 
1 27 
1 28 
000 
1 54 
125 
8 00 
1 
1 
1 
1 
1 
1 
1 



5 
1 
5 
2 



1 
167 7 
158 
000 3 
128 


12 5 
150 
8 
8 

8 2 
2 
1581 
128 5 
128 5 
00 2 



1 
2 
2 
3 
7 
1 
1 
7 



-xi- 

159 
154 
155 
800 
157 



2 
9 
7 
2 
8 J 



13* MULTIPLICATION 

firs, mpy 

V- Get B 



Of5 . NGMPY 
Get -B 



Common steps 



Get A 



If A i-O, C-0 directly. 
A«"Ai,a. Split up and store 
B*Bl,b. Split up and store. 



Calculate and store 
49 - (a+b) 

Go to error stop if exp. of 
prod, would be out of range 

c{. « At. • Bi 



Common with DIV . 

If C{ > 10, C 7 - CJ/10. 



If Cj < 
correct 



10, prepare to 
exponent 



(If Ci-10 (due to 
.rounding), go to 



carry in 
correct. 



Calculate 
attach it 



c • a+b-49-corr. j 
to C], 



-odi- 



13 54 1595 69 5555 133 8 



Error atop, for MPY and DIV. 



1 
1 

1 
1 
1 
1 

1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



3 
3 
3 

3 

3 
3 

3 
3 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 



14 
14 



14 
14 
14 

14 
14 
14 



15 
15 
15 
15 



61 
62 



1 
2 
3 
4 



60 
60 
08 
68 
80 
81 



1 338 
1316 



61 
60 
60 
49 



65 

60 



157 8") 

15 7 8 

16 5 5V Constants 



00 lJ 



8 002 

8 3 




14 4 5 Go to store after error stop, 

15 3 9 

U. DIVISION . 
2 Oj^, DIV 
9) 

Get B 

B«0 tested below 



Get A 



A«*Ai,a. Split up and store 
A,|«0 tested below 



B»Bj »b. Split up and store 



Store shifted A], 



Calculate and store b-a*»50, 

lGo to stop (in MFY) if exp. of 
/quotient would be out of range, 

5l Cj-Ai/Bi; go to exp. corr. 



(in MPT) and storing routines • 



3 3 1610 60 
3 4 1659 60 
3 5 1654 50 



1673 
160 6 



- Constants 



36 1673 
3 7 1527 
3 8 1374 



1350 
1169 
8002 
1137 



4 5 
45 
69 



30 
1 5 
67 
60 



162 6 
168 
3 33 3 



000 3 
10 22 
[ B ] 
8 00 2 



137 A 
144 5 
133 8 



Go to stop if BO. 

Set C-0 directly if Aj-0 

Error stop for B-0. 

15. ABSOLUTE VALUE . 
116 9) 02^2^0, ABS 
8 2 

1 1 3 7 \ Get ISI - C, go to store it, 
144 5j 



-adii- 



15 



5 1022 67 



113 7 Constant 



16 
16 



16 
16 



16 
16 
16 
16 
16 
16 



1 

2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 



2 7 
2 8 



2 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
4 
4 
4 
4 
4 
4 

4 
4 
4 
4 
5 
5 



1300 
1214 
80 02 
1206 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 



39 

48 
29 
51 

3 5 

45 
52 
09 
17 

2 7 
38 
1 6 
25 
82 
36 
91 

4 9 
1 
28 
44 

3 7 



2 22 
20 8 
215 

2 

3 

9 

4 

1 

2 

4 

1 

2 



57 
2 6 



21 9 
2 34 
240 

55 
419 

1 50 



20 
60 



60 
1 9 



60 



1670 



000 
121 
B 
67 
4 44 
8 00 
000 
122 
127 
800 

00 
125 
123 
125 

1 24 
000 
125 
8 
127 
8 01 
1 29 
127 
00 
123 
8 
129 
10 5 
8 02 
8 00 1 
00 
105 
127 


141 

1 27 
8 
124 
115 
127 



4 

2 
2 
3 
8 
3 
1 
5 
4 
5 

1 
5 
2 



8 2 
0002 
1_2 7 8 
1270 
2 



12 14 
800 2 
120 6 



16. SQUARE ROOT 
O y^OO, SCffi T 

Get B 

Stop if B < 0. 

If B*0, C«0 directly. 
B*Bi,b. Split up and store 

x *1+.22B1 



Calc. and store Bj/xjj. 

If [Bt/Xjj - XjJ < 10~ 8 , go to end, 



x^ * 1/2 (Bi/xn + x n ) 



i 



) 



Calc. 1/2 (b+50). 

If b is even, go to end. 

If b is odd, Ci - Xjj vlo. 

Equip with exp, , go to store 
If b is even, Ci*Xjj. 



> Constants 



03 1622 



-xiv- 



17. EXPONENTIAL 
17 1 1301 6 9 1514 1467 Og-301 . EXP E 

17 2 1351 69 1464 1467 Q g »351 . EXP 10 



17 


3 


14 6 7 


30 


00 3 


142 5^ 


17 


4 


1425 


24 


148 5 


148 8 


17 


5 


1488 


15 


1441 


800 2 


17 


6 


80 02 


60 


[ B ] 


1553'} 


17 


7 


15 53 


30 


0002 


1459 


17 


8 


14 59 


21 


126 4 


141 7 J 


17 


9 


1417 


46 


142 


142 1") 


17 


10 


1420 


6 9 


142 4 


147 7 


17 


1 1 


1421 


69 


152 4 


147 7 


17 


1 2 


14 77 


24 


1480 


163 3. 


17 


1 3 


1633 


67 


8 02 


154 li 


17 


1 4 


15 41 


16 


144 4 


155 


17 


1 5 


1550 


35 


001 


1458 


17 


1 6 


14 58 


44 


1462 


1512 


17 


17 


1 512 


30 


00 5 


1 52 5J 


17 


18 


1525 


4 6 


147 8 


152 9") 


17 


1 9 


1478 


69 


1481 


14851 


17 


20 


1529 


69 


148 2 


148 51 


17 


21 


.14 85 


22 


149 5 


155 l 1 


17 


2 2 


1485 


22 


149 5 


150 2. 


17 


2 3 


JL O O JL 


60 


150 4 


1509 


17 


2 4 


1509 


19 


126 4 


163 1 ' 


17 


2 5 


1631 


3 5 


002 


143 7 


17 


2 6 


1437 


67 


8 00 3 


149 5. 


17 


2 7 


1502 


67 


126 4 


166 9) 


17 


2 8 


1669 


35 


2 


149 5j 


17 


2 9 


1495 


[ 1 


[ |b-49l] 
12 7 


1465 


17 


3 


1465 


20 


1474 


17 


31 


1474 


21 


12 7 8 


168 1 


17 


3 2 


16 8 1 


1 1 


158 7 


149 1 


17 


3 3 


1491 


4 6 


159 4 


139 3 


17 


3 4 


1594 


60 


154 8 


151 6} 


17 


3 5 


1516 


1 9 


127 


169 5 


17 


3 6 


169 5 


60 


8 00 3 


165 3 


17 


3 7 


16 53 


10 


150 6 


146 6 


17 


3 8 


1466 


1 9 


1270 


169 4 


17 


3 9 


1694 


60 


8 00 3 


12 5 4 


17 


4 


12 54 


10 


165 8 


1416 


17 


41 


1416 


1 9 


1270 


144 8 


17 


4 2 


1448 


60 


8 00 3 


160 8 


17 


4 3 


1608 


10 


1612 


1667 


17 


4 4 


16 67 


1 9 


1270 


159 8 


17 


4 5 


1598 


60 


8 00 3 


15 5 8 


17 


4 6 


15 58 


10 


1562 


1617 


17 


47 


1617 


1 9 


1 270 


155 2 


17 


4 8 


15 52 


60 


8 00 3 


160 9 


17 


4 9 


1609 


10 


1513 


156 7 


17 


50 


1567 


19 


1270 


140 2 


17 


5 1 


1402 


60 


8 00 3 


155 9 


17 


5 2 


15 59 


10 


146 3 


1517 


17 


5 3 


1517 


1 9 


1270 


160 2 


17 


5 4 


1602 


60 


8 00 3 


166 2 


17 


5 5 


1662 


10 


166 5 


152 
148 OJ 


17 


5 6 


15 20 


19 


8 00 3 



Common steps * 

Set an instr, for EXP E or EXP 10. 

Get B~Bi t b. Split up 
and store. 



Set an instr. distinguishing 
between pos. and neg. exp* 



Form b-49. Go to special 
routines if |b-49| £ 10. 



Set instr. for left or right 
shift of |b~49|. 
Go to EXP E routine 
Go to EXP 10 routine 

For eB-lO 1115 , calc. |mBi|. 

For 1QB, get |Bi|. 

Shift |I|«n+\. (0 ^ \ < 1). 
Store \, 
Store n. 

If |1| > 50, go to stop* 



10^ - [l+a 1 \+.*.+a7\ 7 ] 2 



-XV- 



17 
17 



17 
17 



17 



7 

8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 



7 7 
78 
7 9 
80 
81 



8 8 
89 





1 
2 
3 
4 
5 



480 
687 
4 98 
66 
25 

83 
80 
3J3 
22 
79 
70 
428 
632 
641 
604 
616 
675 
431 
462 
615 

514 
464 
441 
524 
4 24 
4 44 
4 81 
482 
504 
8 7 
4 8 
06 
58 
12 
62 
1 
6 
6 
3 



21 
6 



60 
46 
60 

22 
22 
60 
65 



02 
06 



96 1393 69 




149 5 

149 5 

8 00 3 

X w u D 



429 4 

009 3 

025 5 
1742 

7 29 5 

5 4 3 9 

627 3 

512 9 



For B > 0, round l6\ 
equip with exponent n+50, 
go to store. 



r\ 



For B < 0, calculate 10 



•A 

Equip 10 with exponent 

49-n (if lO"*- < 1) or 

50 (if 10**^ «l)j go to store. 



Go to stop if b > 59, 
set C«l if b £ 39. 



155 r> 

150 2 
155 3 
1687 
143 8 



14 
14 
48 




5 
5 
9 


2 64 3 
4 918 



y Constants 



11 
17 
35 



08 8 4 
27 7 6 

00 5 



6666 1338 Error stop for exp. , logo & trig. 

16. LOGARITHM. 



18 


1 


1302 


69 


1807 


1818 


? «302. LOG E 


18 


2 


1353 


69 


1815 


18 18 


2 -352. LOG 10 


18 


3 


18 18 


24 


1821 


182 4 


Common steps. 


18 


4 


18 24 


30 


000 3 


183 3' 




18 


5 


18 33 


15 


1836 


800 2 




18 


6 


8 02 


6 


[ B ] 


183 1 


- Get B*Bi f b. Split up, 


18 


7 


18 31 


30 


002 


183 7 


store exponent 


18 


8 


18 37 


20 


1291 


184 4. 




18 


9 


18 44 


44 


184 7 


139 3" 


„ Go to stop (in EXP) i 


18 


10 


18 47 


46 


139 3 


17 5 1. 




18 


11 


17 51 


60 


8 00 3 


180 9"] 
1817! 




18 


1 2 


1809 


11 


1812 





-xvi- 



1 8 




8 29 

8 3 8 

8 45 

8 3 

8 81 

8 39 

8 26 

8 88 

8 95 

8 22 

8 8 2 

8 40 

8 10 

8 19 

8 80 

8 35 

8 46 

8 4 9 

8 62 

8 68 



60 
11 

30 



28 
29 
84 
80 
00 
84 
28 
29 



128 5 

000 8 

8 1 

00 4 

183 5 
00 4 

184 9 

186 5 

127 2 



If Bj < e, set u * Bj/Zc 



If Bi > e, set u « Bi/e 1 ^ 5 , 
v«1.65 



u-1 
u+1 



log it£ + 
6e 1-t 



t[2 



|t 2 



3 



t* 



ix 



t 10 ]+v. 



End of common steps 



For LOG E, store L, 
calculate Mb 1# (bj-b-50) 



5 1838) 

1 184 5 
8 180 3 

6 18 8 2 
3 18 3 9 

2 18 2 6 

5 1 8 8 8 I For LOG 10, store mL, 

1 1 8 9 5[ calculate b,. 

8481822 



0001 188 2J 



8 4 0' 


810 


819 


8 8 


89 3 


84 6 


44 5 


86 2 


8 6 8 


87 5 



Common steps 



jGo to stop on loss of more than 
(.two digits of accuracy 
JGo to store machine zero 
I if G ! has seven zeros. 



-XVI 1- 



18 
18 

18 



18 
18 
18 
1 8 

18 

18 
18 
18 



19 
19 
19 

19 
19 
19 



6 9 
70 

71 
72 

7 3 
74 
7 5 
7 6 
7 7 

7 8 
7 9 
81 
82 



84 
8 5 
86 
87 
8 8 

8 9 
80 

9 



1 
2 
3 
4 
5 
6 
7 
8 
9 



00 
01 
02 
3 
04 



1 
2 
3 

4 
5 

6 

7 
8 
9 

1 
2 
3 
4 
5 



8 75 
834 

8 43 



65 
31 



18 4 



811 
8 78 
8 08 
8 79 
8 87 



17 50 

18 54 
1874 
1872 



8 
81 
83 
8 1 
82 
86 
86 
81 
8 3 
8 8 
8 6 
83 
84 
8 
8 4 



7 

5f 

6 

2 

8 

7 

3 

6 



6 

9 

2 

8 

6 

2 

6 

6 

5 

3 

2 

6 

5 



1303 
1614 
1721 

1304 
1664 
1 622 



157 
153 
154 
14 5 
80 
157 
1582 
1689 
12 67 



9 
6 
5 
5 
2 
4 



4 4 
30 
1 5 
16 



83 1893 01 



1 
1 
6 

6 

4 
1 
1 
1 
2 

5 

2 

4 

1 

2 

2 

40 

66 



6 9 
24 
69 

69 
24 
69 



800 3 
002 
002 
180 8 
27 2 
832 
27 2 
832 



8 2 

185 4 
0001 
189 6 
1896 



1274 
127 4 

2 718 
6530 

1872 
5 000 
204 9 



25 
29 
81 
22 
71 



6 66 6 



125 6 

126 7 
1526 

1 707 
1267 
157 6 



128 3 
8 002 
000 3 
1561 

[ B 
127 
00 2 
1252 



83 
84 
75 
81 
87 
8 8 
87 
88 
44 



L+Mb, for LOG E 
mL+b, for LOG 10 



Normalize and round. 




180 4 
186 3 
180 4 
180 4. 

183 5 



82 
88 
83 
81 
5 9 



9] 

1 

1 

8 

7 



In case a carry occurs 
in rounding, correct.* 



Stop on loss of two digits, 
return to log. progr. on start. 



Y 



1874 
908 6 

005 2 

509 3 
4819 
1818 
222 2 
285 7 

666 6 
000 1 



1614 
172 1 

157 9 

1664 
162 2 
1579 

153 6 

154 5 
14 5 51 
800 2 
1574 

158 2 
168 9 
12 6 7 



• Card out of place. 



Constants 



[ ] [ ] 



19.. SINE AND COSINE . 
02*303 . SIN R 



0g«304 . COS R 



Common steps . 

(this step needed in degree prog.) 



Get B»Bi ,b. Split 
up and store. 



acviii- 



19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



19 

19 

19 

19 

19 

19 

19 

19 

1 

1 

1 

1 

1 

1 



1 6 

17 

1 8 

19 

20 

21 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

7 

7 

7 

7 



2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



75 
38 
95 
98 
49 
08 
66 



67 
5 6 



27 

82 
88 
82 
38 
46 




5 
5 
6 
8 
3 
6 
3 
6 
1 
4 
3 



1 

7 

1 

2 

1 

9 

1 

2 

7 

4 
70 
71 
283 
742 

47 

05 

6 

6 

7 

7 



4 

4 

3 

2 

2 



88 
44 



4 

1 

6 

3 

8 

9 

68 

28 
290 
371 
280 



5 
2 
8 

7 
6 



60 
68 
15 



3 

4 

3 

4 

6 

2 

6 

2 

1 

4 



6 

1 

3 



5 

4 



6 

9 

2 

9 

2 

5 

6 

1 



9 

5 



4 6 
10 
61 
10 
10 

y 



4 

1 

4 

1 

4 

1 

1 

6 

3 

1 

4 

2 

6 

1 

2 

6 

1 

6 

1 

1 

6 

1 

1 

6 

1 



8 002 

1691 

129 8 

56 9 

252 

711 

56 9 

27 8 

27 2 

4 4 

00 

27 

00 

68 

63 

29 

53 

29 

170 

137 



1 
1 
1 
1 
1 
1 


1 

1 
1 
1 
1 
1 



4 

1 
6 
5 
2 
5 
3 
5 
3 
2 
5 



1618 
1252 
00 2 



y fe> b i ] 



6 

6 

6 

1 
1 

1 
4 


9 
1 

9 


1 
9 


9 

1 



1 5 3 8 "\ Tests in sine routine : 

1 2 9 5 I 

17 4 9V Test 47-b. 

1674J 

17 8' 

166 6 

167 4 
14 4 5 
162 7 
160 1 



3 

3 

3 

6 

1 

5 

1 

6 

1 

6 

1 

1 

3 

9 

1 

3 

5 

1 

1 



3 


3 
4 



147 6 

188 6 

8 00 

8 00 

8 00 

188 

8 00 

170 

8 00 

126 

8 00 

167 

8 00 

8 00 

8 00 

000 

8 00 

1 38 

128 

8 

8 00 

125 

8 0^0 3 

1572 

8 00 

149 

125 

8 

153 

1250 

8 00 3 

168 4 



If b-47, test |B|-.0025 

If |B| < .0025, go to store C-B 



If |b| > o0025, test 49-b, 



J. W ^ J. I J.J. |U| £ ON~<-y, V~W W -,, ~ 

1 7 9 J go to special routines if 
|49-b|^10. 




For |49~b| < 10, set 
shift instruction. 



1 7 1 0*) If b > 52, stop to indicate 
17 15V loss of accuracy; continue 
5 J on start 

6"] 

Compute the fractional 
V part F of B/2rr, then 

F ^ 0, G«F+1 if F < 0. 




for SIN, 1+4G for COS. 

_ a 1 'a 
17 6} 

12 6 21 if l<a£ 2, a'"2~a 
174 8 J 

12 7 lllf 2< a< 3, 
163 4J 

1540 If 3 < a< 5, 
15 4 0] 

17 4 8 I These angle reductions 
17 2 (carried out, 
174 8 J 

17 2 4llf a'-l, s«t S-l. 
12 8 8] 
17 4 4^ 



a' 



a 1 



*2«a 
"a-4. 



17 
17 
16 
17 
13 
12 
13 
17 
12 
13 
12 
13 



>S«0 1 a l+ 03a l3 +..o*e 9 a 



1 9 



8 
8 9 



-odbc- 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



1 
1 
1 
1 
1 

1 
1 
1 
1 



80 



■8 8 

8 9 

9 
91 
9 2 
9 3 



00 
01 
02 
3 
4 
5 

6 
07 
08 
09 
10 

1 1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



3 8 



7 

3 

9 

2 

4 

7 

7 

8 

2 

8 

9 

79 

37 

63 

25 

87 
2 86 
276 
380 



2 56 
5 26 



7 

6 

9 

1 

3 

3 



02 
18 



19 112 1353 



60 
36 



20 
69 
68 
11 



1 5 



07 
64 
15 



125 
00 
8 
139 
128 
8 00 
00 
138 
000 
00 
000 
00 
8 00 



002 
28 6 
69 
69 
39 3 
86 5 



1272 
188 6 
8 002 
188 6 



7 500 



3 
9 15 4 

15 

61 

68 

96 

07 




4 
9 
5 
7 



49 

69 138 4 



1275 
173 4 
139 4 
172 6 
64 



1 

16 7 2J 

12791 

128 

172 

137 



4 

2 

9 
139 5 

13 7 91 Round and normalize, 
1 7 3 7 /go to store 

3/ 



126 
172 
12 



Go to stop if b ;> 59, 
set a«0 if b £ 39. 



87 
14 4*5 
144 5J 
1380 
128 3 

157 5' 
174 2 
1627 
174 2 
157 4 



171 6 

1-7 1 6 V Constants 



943 
484 2 
765 6 
679 3 
7111 
6318 
005 



17 3 9 Og-353 . SIN D 



19 113 1354 6 9 1385 17 39 2 *354 . COS D 



7 39 
50 5 
585 
698 



2 4 
69 
24 
20 



29 3 
261 
5 3 9 
28 9 



118 1555 60 1167 

119 1384 21 

120 1385 21 

121 1261 35 0002 

122 1167 17 4 532 



150 5) 
158.5 
169 8 J 
15551 
157 8 J 

130 3) 
130 4 
129 3 
934 8 



Common steps 

Set an instr« at the end of 

MPY to get out to SIN R or COS R. 

Go to MPT with B in place of A 

and 2nr/360 in place of B. 



r Constants 



20 



20. ARC TANGENT . 
1 1305 69 1758 1761 0^305 . ARC TAN R 



-XX- 



20 2 13 55 69 1908 1761 Op»355 . ARC TAN D 

Common steps 



Get B-Bi ,bj store Bj • 



Prepare to use the formula 
arctan (-x) ■ -arctan x 
at the end if B < 0. 



20 


3 


1761 


24 


126 4 


1917 


20 


4 


1917 


30 


000 3 


192 5) 


20 


5 


1925 


1 5 


192 8 


800 2 


20 


6 


80 02 


60 


1297 


194 1V 


20 


7 


1941 


21 


190 


20 


8 


1900 


3 


00 02 


17 5 7 


20 


9 


17 57 


21 


126 5 


17 6 8J 


20 


10 


1768 


46 


1771 


177 2) 


20 


IX 


17 7 1 


69 


177 4 


17781 


20 


1 2 


17 7 2 


69 


177 5 


177 8 


20 


1 3 


17 78 


24 


128 3 


17 8 6) 


20 


1 4 


1786 


68 


8 002 


179 5) 


20 


1 5 


1795 


1 5 


156 8 


1773 , 


20 


1 6 


1773 


35 


001 


177 9{ 


20 


17 


1779 


44 


178 3 


178 4 


20 


1 8 


1783 


4 6 


193 7 


193 0J 


20 


19 


17 8 4 


30 


00 5 


189 7) 


20 


20 


18 97 


69 


1901 


175 4V 


20 


21 


17 54 


22 


126 


1763J 


20 


2 2 


1763 


46 


176 6 


1918 


20 


2 3 


1918 


16 


1921 


192 6) 


20 


2 4 


1926 


4 6 


192 9 


193 Of 


20 


2 5 


1930 


60 


129 7 


12 6 4J 


20 


2 6 


1929 


65 


126 5 


17 6 91 
1777J 


20 


2 7 


17 69 


69 


1922 


20 


2 8 


17 66 


60 


1919 


192 3) 


20 


2 9 


1923 


64 


126 5 


177 oy 


20 


30 


1770 


69 


192 4 


177 7J 


20 


31 


17 77 


24 


128 5 


17 8 8) 


20 


3 2 


17 88 


3 5 


0001 


12601 


20 


3 3 


1260 


30 


[49-b] 
8 002 


17 8 1J 


20 


3 4 


17 81 


6 1 


17 8 9) 


20 


3 5 


17 89 


20 


129 3 


179 6 


20 


3 6 


17 96 


1 5 


179 9 


175 2) 


20 


3 7 


17 52 


4 4 


175 5 


175 6) 


20 


3 8 


17 55 


60 


175 9 


191 3) 


20 


3 9 


19 13 


1 9 


129 3 


192 7 


20 


4 


19 27 


3 5 


0001 


193 3 


20 


41 


19 33 


10 


166 5 


192 


20 


4 2 


1920 


21 


127 7 


1 7 8 ' 


20 


4 3 


17 80 


60 


129 3 


179 7 


20 


4 4 


1797 


11 


175 9 


1764 


20 


45 


1764 


64 


,127 7 


190 2 


20 


4 6 


19 02 


6 9 


190 5 


190 9. 


20 


47 


17 56 


65 


129 3 


114 7) 


20 


4 8 


1147 


3 5 


0001 


175 3! 


20 


4 9 


17 53 


6 9 


19 0.6 


190 9 1 


20 


50 


1909 


2 4 


1267 


1 57 OJ 


20 


51 


1570 


20 


127 8 


193 1A 


20 


5 2 


1931 


60 


8 02 


193 9 


20 


53 


19 39 


19 


8 01 


194 5 


20 


5 4 


19 45 


21 


1250 


190 3 


20 


5 5 


1903 


61 


8 00 3 


1911 


20 


5 6 


1911 


1 9 


177 6 


194 


20 


57 


19 40 


60 


8 3 


114 9 


20 


5 8 


1149 


10 


185 6 


1762 



Test 49-b; go to special 
routines if |49-b| > 10 



Prepare for shift of 49-b. 

If b > 49 » go to arctan 1/x. 
If b <J 46, set arctan x • x, 
go to end. 

If 46 < b £ 49, go to 
series evaluation 

If b > 49, calculate l/Bj . 



Shift to get x*B or x«*l/B 
into fixed decimal form 

If jx| > .28, go to calc. 



x 



1 + 



-y 



x y 



; y*.6 



1+ |x|y* 



z ■ |xj . Prepare to 
add or arctan y. 



-xxi- 



30 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 



20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 



8 



8 8 
8 9 
90 



1 
2 
3 
4 
5 
6 
7 
8 
9 







8 
09 



7 
4 
2 
2 
1 



62 

9 
98 

14 

3 8 

48 

60 
7 87 
146 





1 

8 

9 

3 2 

67 
28 5 
943 
910 
934 
9 4-2 

1 5 

9 

1 

3 

4 

6 

6 
2 8 
2 64 
9 36 
937 

758 
908 
928 
74 
75 
68 

2 



3 
6 
5 
4 
5 
4 
3 



7 
7 
5 
9 
9 

9 2 
91 
92 
7 9 
7 5 
66 
T9 
190 
7 8 
7 3 
5 
7 
5 
7 



80 
01 
18 
9 4 



y 



65 



65 



801 



20 

5 

1 5 

57 
50 




128 3 

8 3 
8 00 3 



00 3 
8 3 
0100 
8 00 3 
2 00 



1731 
158 

4 04 1 

7 07 9 

2 95 7 



arc tan z * 
,2 



10 



n It . it it n 

(coefficients from LOG) 



(correction at the end) 

arctan |x|««arctan z+fO or arctan y] 

If b>49, arctan |x| «tt/2- arctan l/|x|. 

If arctan Ib| ^ 1, round 
in 7th decimal place, if 
< 1, round in 8th. 

Normalize arc tan |b| 



C » arc tan I in radians 
For degrees, go to MPY 

by 360/2tr. 

For b > 59, arc tan |b| * n/2. 



144 5 

193 6 

194 1 
126 4 
1264 

178 1 

1910 

1 9 4 3 > Constants 



12 8 5 

12 8 5 

950 

632 7 

005 

805 1 



20 114 1143 19 1801 1267 Correction in arc tan z. 
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